x-—iy 
z 


$65y) tiw sy) —xtiy HR? z 
Xx ty 





atau 


2 


R 
(m3) 


2 
Xx 


R 
van 7 :) 
Xx ty 


Persamaan atur untuk garis arus dengan fungsi arus bernilai nol adalah 


R 
1— -G 
s( 35) 


Solusi dari persamaan diatas adalah 








y-0 
yaitu sumbu-x dan juga 
x in y 2 R 


yaitu suatu lingkaran berpusat di (0,0) dengan jari-jari R. Karena garis arus dapat diganti 
oleh permukaan benda padat yang kedap fluida, maka fungsi potensial kompleks (5.66) 
adalah solusi dari masalah aliran potensial disekeliling suatu silinder bundar dengan 
penampang lintang berbentuk lingkaran. Masalah aliran potensial di sekeliling satu atau 
banyak lingkaran dapat dicoba diselesaikan dengan bantuan Teorema Lingkaran Milne- 
Thomson (1966), yaitu 

Ditentukan bahwa ada aliran potensial, tak-mampat dan 2-dimensional dibidang-z. 
Ditentukan pula bahwa medan aliran tidak dibatasi oleh benda padat kokoh dan fungsi 
potensial kompleks aliran adalah f(z). Semua singgularitas yang ada di medan aliran berada 
di titik-titik dengan jarak yang lebih besar daripada R diukur dari titik asal sistem 
koordinat. Bila sebuah silinder bundar dengan penampang lintang berbentuk lingkaran 
berjari-jari R, yaitu lzl - R, diletakkan dalam medan aliran f(z) maka potensial kompleks 


aliran berubah menjadi F(z), dimana 


F(z)- tas) 
zZ 


Sebagai contoh, bila dalam suatu aliran angin seragam dengan fungsi potensial kompleks 
f(z2)—z.exp(—ia) 
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kita letakkan lingkaran |z| - R, maka fungsi potensial kompleks f(z) berubah menjadi F(z) 


dimana 


F(z2)- sean et 
z 


Khususnya apabila 4-0, maka 


2 
Ftej ol 
z 


Perhatikan bahwa rumus untuk Ffz) itu tepat sama seperti persamaan (5.66). Masalah 
aliran di sekeliling lingkaran akan dibahas secara rinci di bab berikutnya. 
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BAB VI 
ALIRAN POTENSIAL DI SEKELILING LINGKARAN 


Tujuan akhir buku ini adalah mempelajari pola aliran potensial di sekeliling bentuk 
aerofoil sembarang. Pengetahuan yang diperoleh kemudian dapat dimanfaatkan dalam 
analisis dan perancangan aerodinamika bentuk-bentuk aerofoil yang diinginkan. Dalam bab 
ini kita akan membahas masalah aliran potensial di sekeliling bentuk lingkaran, karena 
masalah ini dapat diselesaikan dengan mudah. Disamping itu aliran ini dapat digunakan 
untuk memperkenalkan konsep syarat batas dan beberapa konsep aerodinamis seperti gaya 
angkat, gaya tahan dan lain sebagainya. Selanjutnya pola aliran di sekitar lingkaran dapat 
ditransformasikan menjadi aliran di sekitar aerofoil dengan bantuan teori pemetaan 
konformal. 

Fungsi potensial kompleks untuk aliran angin seragam yang menyapu lingkaran 
(el-R di bidang-z adalah 

O(z)-24#-R'/z (6.1) 
Persamaan di atas dapat diuraikan menjadi persamaan untuk fungsi potensial kecepatan 
dan fungsi arus sebagai berikut 


d(x,y) —x4-R7/ (4) (6.2) 

v(x,y) 5y(-R/(x 4y”)) (6.3) 
Persamaan atur untuk garis-garis arus aliran ini dapat diperoleh dari persamaan (6.3) 
sebagai berikut 

Y WAY # (PR )y- y,x7 50 (6.4) 


dimana w, dan R adalah konstanta tertentu. Persamaan (6.4) dapat ditulis dalam sistem 
koordinat polar menjadi 
r'sin0-—w,r—R”sin0 —0 (6.5) 
PX 4y 
tan0 — y/x 
Koordinat titik-titik pada garis arus dengan harga fungsi arus w,, dapat dihitung sebagai 
berikut. Apabila titik P dengan koordinat (XP,YP) berada pada garis arus W,, maka kita 
dapat menentukan harga XP, dan YP adalah solusi dari persamaan berikut 
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yY Ay 4By4-C-—0 (6.6) 


dimana 


A--—y, 

B-XP'-R? 

C--—y,XP? 
Sebaliknya, apabila harga YP ditentukan maka XP adalah solusi dari persamaan berikut 

xx —D/E (6.7) 
dimana 

D- (R' 4 y,YP- YP?)YP 

E-YP-y, 


Dalam koordinat polar, koordinat titik P adalah (RP, OP) dan apabila 0P ditentukan 
maka RP adalah solusi dari persamaan 


r—A#vA?'4R? (6.8) 


dimana 
A —y, / 2sinOP 
Sebaliknya apabila harga RP ditentukan maka OP adalah solusi dari persamaan berikut. 


RP 


La Du (6.9) 
Khusus untuk y, — 0, persamaan (6.4) dapat ditulis ulang sebagai berikut 

(24 -R?)-0 (6.10) 
dengan solusi 

y-0 (6.11) 

4 YR? (6.12) 


Solusi di atas menunjukkan bahwa garis arus w,, — 0 adalah garis sumbu x dan lingkaran 
lzl - R. Pola aliran yang diberikan oleh fungsi arus (6.3) dapat dilihat pada gambar (6.1). 


Kecepatan fluida di suatu titik dalam medan aliran dapat diperoleh dari rumus berikut 


O'(2)-u-iv—1-(R/z)7 (6.13) 
atau 

u(x,y)-1-— (2 —y).R?/ (2 4y) (6.14) 

V(xry) 5 —2xyR/ (2 4y (6.15) 
Komponen kecepatan dalam sistem koordinat polar adalah 

V.(r,0) — (A-R?/r”)cos9 (6.16) 
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V,(r,0) —-—(1#R?/r”)sin0 (6.17) 
Pada permukaan lingkaran, yaitu untuk 2 4#y'-R?, harga komponen-komponen 


kecepatan adalah 


u(x,VR? —x)- 20-— X/R)) 
VX, VR? —x)- —2x/R'—x /R? 


V.(R,0) -—0 
V(R,0) — 2 
Hasil-hasil tersebut di atas dijabarkan sebagai berikut, pada permukaan silinder bundar 

(lingkaran) aliran bergerak arah tangensial atau sesuai arah garis singgung pada lingkaran. 

Kecepatan tersebut berharga negatif untuk belahan atas lingkaran, yaitu 0 « 0 «mx, karena 


| untuk -RXxSR (6.18) 


untuk 0x0x27 (6.19) 


arah kebalikan putaran jarum jam didefinisikan sebagai positif dalam perjanjian matematis. 
Pada belahan atas lingkaran fluida mengalir dari kiri ke kanan sesuai arah putaran jarum 
jam dan kecepatan bernilai negatif. Sebaliknya pada belahan bawah lingkaran, yaitu 
mXO« 2x, kecepatan tangensial bernilai positif karena aliran fluida dari kiri ke kanan itu 
sesuai arah kebalikan putaran jarum jam. 

Komponen kecepatan radial selalu berharga nol dan ini sesuai dengan syarat batas 
bahwa fluida tidak boleh masuk menembus permukaan silinder yang dianggap kedap 
fluida. Pada titik (-R,0) kecepatan aliran berharga nol. Pada saat partikel fluida menabrak 
titik tersebut dia terhenti dan kecepatannya menjadi nol. Karena partikel fluida di titik 
tersebut menjadi stagnan atau diam tak bergerak, maka titik (-R,0) disebut titik stagnasi. 

Aliran juga dapat dibayangkan sebagai fluida yang diam dan silinder bergerak dari 
kanan ke kiri. Titik (-R,0) dengan demikian berada di bagian depan silinder dan disebut 
titik stagnasi depan. Titik itu juga dikenal sebagai titik hidung karena berada di ujung 
depan. Kecepatan fluida meningkat dari titik hidung sepanjang permukaan silinder dan 
mencapai harga maksimalnya di titik pundak atas, (O,R), dan di titik pundak bawah, (0,- 
R). Harga maksimal tersebut adalah 2, yaitu dua kali kecepatan arus bebas atau dua kali 
kecepatan gerak silinder di udara yang diam. Berawal dari titik pundak, kecepatan aliran 
secara berangsur-angsur diperlambat sampai akhirnya mencapai harga minimal, yaitu nol, 
di titik (R,0) yang disebut titik stagnasi belakang atau titik ekor. Fluida yang meninggalkan 
titik ekor kemudian dipercepat gerakannya sampai akhirnya mencapai harga arus bebas 
yang tidak terpengaruhi oleh adanya silinder bundar. Distribusi kecepatan di sepanjang 
sumbu-x dan lingkaran, yaitu garis arus w-0, yang dijabarkan di atas diberikan dalam 
bentuk grafis pada gambar (6.2). 


69 


Garis arus dimana titik-titik stagnasi berada, untuk kasus ini adalah garis arus w — 0, 
disebut garis arus stagnasi. Karena garis arus tersebut membelah aliran menjadi belahan 
atas dan belahan bawah, maka garis tersebut juga dikenal sebagai garis arus pembelah. 
Persamaan Bernoulli menyatakan bahwa tekanan total di sepanjang garis arus itu bernilai 
tetap. Karena di kejauhan tak-berhingga dari silinder tekanan statik bernilai P,, dan 
kecepatannya adalah V,,, maka persamaan Bernoulli dapat ditulis sebagai 


1 1 2 
Pt —pV'—- Pt —pV, 
2 & 2 & 


Koefisien tekanan adalah tekanan statik tanpa dimensi, dengan simbol C, dan dirumuskan 


sebagai berikut 


CP - Bo) /TPVe SI (VI Va)” 


atau dalam notasi kecepatan tanpa dimensi (lihat bab V) 


C,-1-V (6.20) 
Distribusi koefisien tekanan di sepanjang garis arus stagnasi diberikan oleh rumus berikut 
C, -1-4sin?'9 untuk 0x0x27 (6.21) 

R' R 
Cc, Be Lan untuk xS-—R dan x2R (6.22) 


Distribusi C, pada sumbu-x dan pada permukaan lingkaran diberikan dalam bentuk grafis 


pada gambar (6.3) dan (6.4). Arah tekanan adalah selalu tegak lurus pada permukaan 
dimana tekanan beraksi. Grafik pada gambar (6.4) dengan jelas menunjukkan kenyataan 
tersebut dengan menggunakan simbol anak panah untuk mewakili tekanan. Tekanan 
negatif digambar sebagai anak panah tegak lurus pada lingkaran dengan arah menjauhi titik 
pusat lingkaran. Sebaliknya tekanan positif digambar sebagai anak panah yang menuju ke 
titik pusat. Konvensi atau perjanjian penggambaran ini sesuai dengan kanyataan bahwa 
dampak dari distribusi tekanan adalah mendorong benda cenderung bergerak ke arah yang 
sesuai dengan arah anak panah. Tekanan negatif cenderung menyedot benda, atau menarik 
benda ke arah kawasan bertekanan negatif. Sebaliknya tekanan positif cenderung 
mendorong benda menjauhi kawasan bertekanan positif tersebut. 

Tekanan adalah gaya per satuan luas, jadi distribusi tekanan dapat dijumlahkan atau 
diintegralkan untuk menghitung gaya resultante atau gaya netto yang beraksi pada silinder 
bundar. Gaya resultante tersebut selalu dapat diuraikan menjadi 2 komponen gaya yang 
saling tegak lurus. Komponen pertama adalah gaya yang beraksi arah tegak lurus pada 
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arah gerak benda atau arah arus bebas. Komponen gaya ini cenderung mengangkat benda 
dan disebut gaya angkat. Komponen gaya yang kedua beraksi ke arah sejajar dengan arah 
arus dan cenderung menghambat atau menahan pergerakan benda, jadi disebut gaya tahan. 
Gaya angkat (lift) dan gaya tahan (drag) yang beraksi pada silinder bundar dapat dihitung 
sebagai berikut. Pertama-tama kita tentukan bahwa panjang silinder adalah 1 satuan 
panjang. Ini berarti bahwa gaya-gaya yang dihitung di bawah ini adalah gaya-gaya per 
satuan panjang silinder. Elemen gaya yang beraksi pada segmen lingkaran R.d@ adalah 
dF—-P.R.do.? (6.23) 

yang dapat diuraikan menjadi komponen gaya angkat, dL, dan komponen gaya tahan, dD, 
sebagai berikut 


dF - dD.£4dL.j (6.24) 
dimana 

dD --P.R.d9.cos@ (6.25) 

dL —-P.R.dO.sin0 (6.26) 


Gaya tahan dan gaya angkat dapat diperoleh dengan mengintegralkan kedua persamaan di 
atas sebagai berikut 


D- -PP.R.cos0d0.x (6.27) 


L-—$P.R.sin0dO.7 (6.28) 
Dari teori kalkulus dasar diketahui bahwa 


PP..R.cos0d0 -0 


9e..R.sinOdO — 0 


dimana P,, dan R adalah konstanta. Ini berarti bahwa persamaan (6.27) dan (6.28) dapat 
ditulis ulang sebagai berikut 


D--$(P- P,).R.cos0d0 


L- be —P,).R.sin0dO 





Mengingat bahwa 
CG: Ta -1-—4sin'0 
—PV., 
2 


11 


maka kita dapat merumuskan koefisien gaya tahan dan koefisien gaya angkat, yaitu C, dan 


C,, sebagai berikut 


Cc, » DM $G,.cosoa0 (6.29) 
1 2 3, 
—PV, -€ 
2 

Teja ——- — 1 fc,.sindao (6.30) 
—PV,”.c 2 


2 
dimana c adalah panjang acuan, yaitu talibusur silinder atau diameter lingkaran, c - 2R. 


Pemasukan rumus untuk C, akhirnya memberikan hasil-hasil berikut 


1 4 
C,- -5pa- 4sin?0).cos0.d0 — 0 


G- — HA Asin? 0).sinO.d0 -9 


dan tentu saja 

D-0 dan L-0 
Seperti dapat dilihat, untuk kasus aliran potensial yang dipelajari ternyata gaya tahan dan 
gaya angkat yang beraksi pada silinder berharga nol. Kenyataan bahwa gaya angkat itu 
berharga nol dapat dimengerti sesuai dengan penalaran wajar, karena distribusi koefisien 
tekanan itu simetri terhadap sumbu-x. karena tekanan di belahan atas silinder itu sama 
dengan tekanan di belahan bawah silinder, maka gaya netto arah sumbu-y, yaitu gaya 
angkat, jelas berharga nol. Kenyataan bahwa gaya tahan juga berharga nol, sebaliknya 
sangat sulit dimengerti karena bertentangan dengan pengamatan fisis yang dapat dilakukan 
sehari-hari. Sebuah silinder bundar, misalnya balok kayu, yang dibenamkan di aliran air 
sungai pasti akan hanyut dibawa arus. Bila kita mencoba mencegah kehanyutan silinder 
tersebut, misalnya dengan mengikatnya dengan tali, kita pasti akan merasakan suatu 
tarikan pada tali yang berasal dari gaya tahan yang beraksi pada balok kayu tersebut. 

Perbedaan antara hasil teori dan pengamatan fisis ini dikenal sebagai Paradoks 
D'Alembert dan menghambat perkembangan ilmu Dinamika Fluida. Masalah ini baru dapat 
diatasi setelah Prandtl (1904) mengemukakan dan menulis tentang teori lapisan batas. 
Menurut Prandtl, walaupun viskositas fluida pada umumnya bernilai sangat kecil namun 
tetap memegang peran yang sangat penting dalam teori pergerakan fluida, khususnya 
dalam suatu lapisan tipis dekat permukaan benda. Dalam kenyataannya semua fluida itu 
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melekat pada permukaan benda, yaitu kecepatan partikel fluida pada permukaan benda 
selalu berharga nol. Sedikit agak jauh dari permukaan benda, melintasi kawasan tipis yang 
disebut lapisan batas, kecepatan fluida mencapai harga hasil prediksi teori aliran potensial. 
Di dalam lapisan batas distribusi kecepatan fluida berubah dari nol pada permukaan 
menjadi kecepatan aliran potensial di luar lapisan batas tersebut. Ini berarti bahwa di dalam 
lapisan batas terdapat suatu distribusi gradien lintang kecepatan yang bernilai sangat besar. 
Hukum viskositas Newton menyatakan bahwa gaya geser yang beraksi pada partikel fluida 
itu sama dengan viskositas dikalikan dengan gradien lintang kecepatan. Walaupun 
viskositas bernilai sangat kecil, tetapi bila gradien lintang kecepatan bernilai sangat besar 
maka gaya geser yang beraksi tidak akan bernilai nol atau kecil, tetapi dapat bernilai cukup 
besar. Inilah sebabnya mengapa gaya tahan, yang merupakan integral dari tegangan geser 
pada permukaan benda, tidak bernilai nol dan dapat bernilai cukup besar. 

Dalam teori aliran potensial permukaan benda dianggap licin dan partikel fluida selip 
atau menggelincir pada permukaan tanpa hambatan tegangan geser. Inilah sebabnya 
mengapa teori aliran potensial memberikan hasil yang bertentangan dengan penalaran 
wajar, yaitu bahwa gaya tahan selalu bernilai nol. Namun demikian ini tidak berarti bahwa 
teori aliran potensial itu tidak ada gunanya karena teori ini dapat dimanfaatkan untuk 
menghitung gaya angkat dengan cukup akurat. Disamping itu, dengan bantuan teori 
lapisan batas, teori aliran potensial dapat juga dimanfaatkan untuk menghitung gaya tahan. 

Sekarang marilah kita lakukan suatu eksperimen otak dalam pikiran kita. Silinder 
bundar yang dipelajari tadi sekarang diputar dengan kecepatan putar @. Karena fluida 
melekat pada permukaan maka kecepatan partikel fluida tentu sama dengan kecepatan 
permukaan yaitu @. R arah tangensial. Gerak berputar silinder tersebut tentu saja juga 
berdampak pada partikel fluida yang berada agak jauh dari permukaan. Dapat diperkirakan 
bahwa dampak tersebut akan menjadi semakin lemah pada titik-titik yang semakin jauh 
dari permukaan. Apabila silinder diputar dengan arah sesuai putaran jarum jam, maka 
kecepatan fluida di belahan bawah akan diperlambat, sedangkan di belahan atas aliran 
menjadi dipercepat. Letak titik-titik stagnasi dapat diperkirakan turun ke bawah. Karena 


kecepatan V meningkat di belahan atas, maka C, -1- V? akan menjadi lebih negatif 
sedangkan di belahan bawah nilai C, akan menjadi kurang negatif. Gaya tarik ke atas oleh 
C, negatif di belahan atas akan menjadi lebih besar sedangkan gaya tarik di belahan bawah 


akan mengecil. Dampak keseluruhannya adalah bahwa gaya tarik ke atas oleh tekanan 
negatif di belahan atas menjadi lebih besar daripada gaya tarik ke bawah oleh tekanan 
negatif di belahan bawah dan ada gaya netto yang cenderung mengangkat silinder. Ini 
berarti bahwa bila silinder diputar arah putaran jarum jam maka suatu gaya angkat positif, 
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atau arah ke atas, akan tercipta. Dipandang dari segi aliran fluida di luar silinder bundar 
sesungguhnya aliran ini mirip sekali dengan aliran tornado yang kemudian diidealisasikan 
sebagai aliran vorteks. Kekuatan vorteks adalah sirkulasi, yang dapat dihitung pada 
permukaan silinder yaitu sebesar 27@R”. Jadi aliran angin seragam yang menyapu sebuah 
silinder berputar, dengan kecepatan putar @, dapat dimodelkan secara matematis sebagai 
kombinasi linear antara angin seragam, doblet danm vorteks. Fungsi potensial kompleks 
untuk aliran ini dalam konvensi (perjanjian) matematis adalah 
D(z)- PAP ja az (6.31) 
z 21 
atau dalam konvensi aerodinamis sebagai berikut 


2 

din) st it Rnr (6.32) 
z 21x 

Untuk menghindari penggunaan harga y yang negatif bila konvensi matematis yang 

digunakan, maka dalam pembahasan di bawah ini kita akan menggunakan konvensi 

aerodinamis atau persamaan (6.32). Perlu diingat bahwa dalam konvensi aerodinamis y 

bernilai positif berarti arah sesuai putaran jarum jam. Fungsi potensial kecepatan dan 


fungsi arus untuk masalah aliran yang dibahas adalah 





#am salin 5. )-a(1) 





KEY 2x x 
2 (6.33) 
R R 
vam rda z 2 ane 13) 
xX ty 41x 
atau dalam koordinat polar 
2 
&(1,0) — ema “5) Rp 
r 2x 
(6.34) 


R: 
Y(r,0) — pain (SJ ir 
r 2x 
Komponen-komponen kecepatan dapat diperoleh dari rumus untuk turunan pertama O, 
yaitu 





2 
R ti 1 


2 


D'(z)-u-iv-1 : 
z 21 z 
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atau 


Heni 
u(x,y)—1-R? 2 rx. ta Pn 
(X4-y)Y 21 x #ty 











: si (6.35) 
2 xy Ko x 
“sy) (24) 217 24y 
Komponen kecepatan dalam koordinat polar adalah 
V, (1,0) — (A-(R/r)).cos0 
(70) (WR /1P)-c05 Wi 
Vp(r,0) 5 —(1-4#-(R/r)).sin0 — YR /27r 


Sekarang marilah kita periksa apakah persamaan (6.32) memang benar-benar mewakili 
fungsi potensial kompleks untuk aliran potensial di sekeliling lingkaran dengan persamaan 
atur x #y” — R?. Pemasukan koordinat (x,y) untuk titik yang berada pada lingkaran ke 
dalam rumus untuk fungsi arus (6.33) memberikan hasil berikut 


w(x,y) —w(R,0)— TEMR (konstanta) (6.37) 
T 


Hasil di atas menunjukkan bahwa lingkaran lzl - R memang adalah sebuah garis arus, jadi 


persamaan (6.32) adalah fungsi potensial kompleks aliran di sekeliling lingkaran tersebut. 
Komponen-komponen kecepatan pada permukaan lingkaran adalah 


utan 32 
untuk x4y -R? (6.38) 

x yR 

ey) Ea Ai. 


Dalam koordinat polar, komponen kecepatan adalah 


V,(R,0)-0 
3 Y untuk 0x0x217 (6.39) 

V, (R,0) — —| 2sin9 - — 

21 
Komponen kecepatan radial ternyata bernilai nol, jadi memenuhi persyaratan bahwa udara 
tidak boleh menembus masuk atau keluar lingkaran. Vektor kecepatan diberikan oleh 
keceaptan tangensial dan bernilai nol di titik stagnasi depan dengan koordinat (R,1 --0,,) 
dan di titik stagnasi belakang dengan koordinat (R,21— 0,,) dimana 
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@, — sin” (y 141) (6.40) 
Garis arus stagnasi tersebut membelah medan aliran menjadi dua bagian yaitu belahan 
medan aliran atas dan belahan medan aliran bawah. Lingkaran juga terbagi menjadi belahan 
atas dengan koordinat 


—0 ,s0xmtO, (6.41) 
dan belahan lingkaran bawah dengan koordinat 
m0, SOx217—0, (6.42) 


Vektor kecepatan di sepanjang belahan atas lingkaran bernilai positif, sesuai dengan 
konvensi aerodinamis yaitu putaran searah putaran jarum jam adalah positif. Di belahan 
bawah kecepatan bernilai negatif karena aliran dari kiri ke kanan memberikan putaran 
kebalikan arah putaran jarum jam. Kecepatan mencapai nilai maksimalnya di titik pundak 
atas atau titik (R,1/ 2), yaitu 
Voaratas “ 2-4#Y/ 21 (6.43) 
dan nilai maksimal dibelahan bawah terjadi di titik pundak bawah dengan koordinat 
(R,31/2), yaitu 

Vinarybawan 5 2—Y / 21 (6.44) 
Harga 0,, atau letak titik-titik stagnasi jelas sangat tergantung pada kekuatan vorteks y . 
Harga y terbesar yang diperbolehkan adalah 

Y mas 5 4T (6.45) 
Untuk kasus ini 0, -1/2 dan kedua titik stagnasi menjadi berhimpit di titik pundak 
bawah, (R,31/2). Seluruh permukaan lingkaran menjadi belahan atas lingkaran, 
sedangkan belahan bawah lingkaran hanya terdiri dari satu titik saja yaitu titik stagnasi 
depan dan belakang yang berhimpit di titik pundak bawah. Bila y bernilai lebih besar 
daripada Y,,..,, maka garis arus stagnasi di sebelah kiri lingkaran akan berpotongan dengan 
garis arus stagnasi di sebelah kanan lingkaran. Ini jelas melanggar pengertian fisis tentang 
garis arus sebagai trajektori pergerakan partikel fluida. Koordinat titik-titik di sepanjang 
suatu garis arus dapat dihitung dengan menyelesaikan persamaan 


y bim? 4y)- va je (X2 “Royal Mint 4-y)- va je -0 
4x 4x 
(6.46) 
dimana y, W, dan R adalah konstanta yang diketahui. Titik P dengan koordinat (XP,YP) 
pada garis arus W,, dapat ditentukan letaknya dengan menetapkan harga x — XP dan YP 
diperoleh sebagai solusi persamaan (6.46) atau dengan menetapkan harga y - YP dan 
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menghitung XP sebagai solusi persamaan (6.46). Persamaan (6.46) adalah suatu contoh 
dari persamaan transendental yang dapat diselesaikan secara numerik walaupun solusi 
analitisnya tidak diketahui. Dalam sistem koordinat polar rumus untuk garis arus diberikan 
oleh persamaan berikut 


R'. YR 
-£ 
r? ) 2 








r.sin@.(1 Inr—y, 50 (6.47) 


T 
yang dapat ditulis ulang sebagai 
2 
ing 5, Ar (6.48) 
2x r 

Apabila titik P dengan koordinat (RP, OP) berada pada garis arus W,, maka RP dapat 
ditetapkan dan harga OP dapat dihitung dengan mudah menggunakan rumus (6.48). 
Bentuk pola aliran dengan harga y — dapat dilihat pada gambar (6.5), diberikan 


sebagai vektor-vektor kecepatan di sepanjang beberapa garis arus tertentu. 
Distribusi kecepatan dan koefisien tekanan di sepanjang permukaan lingkaran 
diberikan oleh rumus-rumus berikut 


V(R,0) - —(2sin0 ——) 
ti untuk 00X (6.49) 
C,(R,0) 1- (2sin0 4-1) 
1 


Koefisien gaya tahan dan koefisien gaya angkat dapat dihitung dengan rumus (6.29) dan 
(6.30) dengan hasil berikut 
Cik — fli-(esino -y 12x) |cosdd0 


C,-0 (6.50) 


C- —li— (asin ty /2n)' JsinOdO 


C,-y (6.51) 
Seperti dapat dilihat, koefisien gaya tahan berharga nol. Hasil ini sudah dapat diduga 
sebelumnya karena dalam teori aliran potensial gaya tahan selalu berharga nol. Tetapi 
perhatikan bahwa koefisien gaya angkat untuk kasus ini ternyata mempunyai suatu harga 
positif tertentu, yaitu dua kali kekuatan vorteks atau sirkulasi tanpa dimensi. Sirkulasi 


1 


mempunyai dimensi kecepatan kali panjang dimana satuan kecepatan adalah V., m/detik 
sedangkan satuan panjang adalah jari-jari lingkaran atau R meter. Gaya angkat dapat 
diperoleh dari rumus (6.30), yaitu 


L- HPV. — pV., R.C, 


atau 

L-— PV, R.y 5 YV,,:PV...R 
dimana c adalah panjang acuan, yaitu R meter. Kuantitas y adalah sirkulasi tanpa dimensi. 
Karena kecepatan acuan adalah V,, m? /detik, sedangkan panjang acuan adalah R meter, 


maka sirkulasi berdimensi dengan satuan m? /detik adalah 


T-yV,.R (m? /detik) 
Rumus untuk gaya angkat dalam variabel berdimensi sekarang dapat ditulis sebagai berikut 
L-pV,Tr N) (6.52) 


Rumus di atas adalah gaya angkat yang beraksi pada silinder bundar per satuan panjang 
silinder. Persamaan (6.52) adalah pernyataan matematis dari teori yang dikenal sebagai 
Teori Kutta-Joukowski, yang menyatakan bahwa gaya angkat per satuan panjang silinder 
itu berbanding lurus dengan sirkulasi di sekeliling silinder. Walaupun tidak dibuktikan di 
sini, tetapi perlu dicatat bahwa teori Kutta-Joukowski di atas itu berlaku untuk silinder 
dengan penampang lintang berbentuk sembarang dan tidak terbatas pada bentuk lingkaran 
saja. Apabila harga p dan V., sudah ditentukan, maka harga gaya angkat, L, hanya 
tergantung pada harga sirkulasi, ', saja. Karena dalam aerodinamika kita berusaha 
mendapatkan gaya angkat yang sebesar-besarnya, maka kita harus mencari bentuk-bentuk 
aerofoil dengan harga sirkulasi yang sebesar mungkin. Bentuk aerofoil yang baik adalah 
yang mampu menghasilkan gaya angkat yang besar dan gaya tahan yang kecil, atau nisbah 
gaya angkat dibagi gaya tahan yang sebesar-besarnya. Gaya tahan muncul karena adanya 
tegangan geser yang beraksi pada permukaan benda dan karena terjadinya pelepasan 
aliran, yang membuat pola aliran menjadi jauh berbeda dari hasil prediksi teori aliran 
potensial. Tegangan geser dan pelepasan aliran itu terjadi sebagai dampak dari kenyataan 
bahwa viskositas fluida itu walaupun bernilai kecil tetapi bukan nol. Untuk mendapatkan 
gambaran sedikit tentang dampak viskositas, di bawah ini kita akan membahas secara 
singkat teori lapisan batas, walaupun hanya secara kualitatif. 

Distribusi kecepatan dan koefisien tekanan pada titik-titik sepanjang permukaan 
lingkaran dapat dihitung dengan bantuan persamaan (6.49). Hasil hitungan untuk kasus 
dimana y -1/2-7 ditunjukkan pada gambar (6.6). Bentuk kurva distribusi C, ternyata 
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seperti mangkok, dan lebih curam untuk belahan atas lingkaran bila dibandingkan dengan 
untuk belahan bawahnya. Pergerakan partikel fluida dipengaruhi oleh 2 faktor yaitu 
kecepatan gerak dan tekanan statik yang beraksi pada partikel. Untuk mendapatkan 
pengertian tentang pergerakan partikel fluida di sepanjang permukaan lingkaran yang 
dihambat oleh tegangan geser viskos, pertama-tama kita bahas dulu suatu masalah lain 
yang ada kemiripannya. Partikel fluida yang bergerak dari titik stagnasi depan ke titik 
stagnasi belakang dapat diibaratkan sebagai pergerakan sebuah kelereng. Kelereng 
tersebut bergerak dari bibir kiri mangkok (titik stagnasi depan) jatuh menggelinding ke 
dasar mangkok (titik pundak lingkaran) dan kemudian naik mendaki sisi kanan mangkok 
sampai akhirnya mencapai bibir kanan mangkok (titik stagnasi belakang). Pada awalnya, 
saat berada di bibir kiri, kelereng berkecepatan nol dan bergerak menuruni dinding 
mangkok karena adanya selisih energi potensial atau gaya tarik bumi. Partikel fluida 
bergerak dari titik stagnasi depan ke arah pundak lingkaran karena adanya selisih tekanan 
statik. Tekanan statik yang besar mendorong partikel fluida bergerak di sepanjang 
permukaan lingkaran menuruni mangkok kurva distribusi tekanan, menuju ke dasar 
mangkok yaitu titik pundak lingkaran dimana tekanan statik bernilai minimum. Kelereng 
bergerak semakin cepat waktu menuruni lereng mangkok karena energi potensialnya 
menjadi semakin kecil. Energi potensial adalah fungsi ketinggian dan energi potensial 
kelereng menjadi minimum saat kelereng berada di dasar mangkok. Dari hukum kekekalan 
energi diketahui bahwa energi tidak dapat diciptakan atau dimusnahkan. Untuk kasus 
pergerakan kelereng, sesungguhnya energi total kelereng itu bernilai tetap walaupun 
kelihatannya kelereng kehilangan energi potensialnya. Energi total kelereng adalah jumlah 
dari energi potensial dan energi kinetik. Energi potensial kelereng, yang kelihatannya 
hilang waktu kelereng menuruni sisi mangkok, sesungguhnya diubah menjadi energi 


kinetik, mv), yaitu kelereng bergerak semakin cepat. Dalam pergerakan partikel fluida, 


energi total partikel itu adalah jumlah energi potensial yang diwakili oleh tekanan statik 
dan energi kinetik. Gerak partikel fluida menjadi semakin cepat karena penurunan tekanan 
statik diimbangi oleh tekanan dinamik (energi kinetik) yang meningkat. Keadaan yang 
dijelaskan adalah idealisasi dari kenyataan sebenarnya. Kelereng yang menggelinding ke 
dasar mangkok sesungguhnya dihambat geraknya oleh gaya gesek. Sebagian dari energi 
total kelereng itu diubah menjadi panas melalui mekanisme gesekan antara kelereng dan 
mangkok. Energi potensial dapat diubah menjadi energi kinetik dan sebaliknya energi 
kinetik dapat diubah menjadi energi potensial. Baik energi potensial atau energi kinetik 
dapat juga dengan mudah diubah menjadi energi panas tetapi energi panas tidak dapat 
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diubah menjadi energi potensial atau kinetik. Energi total kelereng sesungguhnya terdiri 
dari 3 jenis energi, yaitu energi potensial, kinetik dan panas. Energi panas itu selalu 
meningkat, jadi jumlah energi potensial dan kinetik dengan sendirinya selalu berkurang. 
Karena adanya disipasi panas, yaitu pengubahan energi jenis lain menjadi panas, maka 
tidak seluruh kehilangan energi potensial diubah menjadi energi kinetik. Kelereng yang 
bergerak dari bibir kiri mangkok akan dapat menggelinding sampai ke bibir kanan 
mangkok seandainya tidak ada disipasi panas. Tetapi karena disipasi panas itu terjadi maka 
kelereng tidak akan mampu mencapai bibir kanan mangkok dan akan berhenti di satu titik 
pada lereng kanan mangkok. Disipasi panas juga terjadi dalam pergerakan partikel fluida 
sebagai dampak dari viskositas. Tegangan geser yang muncul di lapisan batas akan 
menentang atau menghambat pergerakan fluida. Dalam kenyataannya kecepatan partikel 
fluida pada permukaan silinder akan menjadi nol sebelum mencapai titik stagnasi belakang. 
Di luar lapisan batas tegangan geser berharga nol, jadi hukum Bernoulli itu berlaku. 
Partikel fluida di luar lapisan batas bergerak cepat dan cenderung menyeret lapisan fluida 
di bawahnya, dalam lapisan batas, yang bergerak lebih pelan. Pergerakan partikel fluida di 
lapisan batas memang jauh lebih rumit daripada pergerakan kelereng dengan gaya gesek. 
Karena adanya hambatan tegangan geser viskos, ada kemungkinan bahwa partikel fluida 
bergerak arah kebalikan dari arah sebelumnya setelah melewati suatu titik tertentu. 
Penyebabnya adalah tekanan yang lebih besar di belakang mendorong partikel fluida 
bergerak maju ke arah depan dimana tekanan statiknya berharga lebih kecil. Karena fluida 
di dekat permukaan bergerak ke kiri sedangkan di luar lapisan batas fluida tetap bergerak 
ke arah kanan, maka dapat dibayangkan adanya suatu aliran berputar atau pusaran di 
kawasan tersebut. Kawasan seperti itu disebut kawasan gelembung pusaran dimana 
partikel-partikel fluida bergerak berputar-putar dilokasi tertentu. Proses terjadinya 
gelembung pusaran dapat dicoba dijelaskan lebih lanjut dengan bantuan sketsa-sketsa pada 
gambar (6.8). Gambar (6.8.2) menunjukkan sebuah partikel fluida berbentuk segi panjang 
dalam suatu aliran potensial dengan arah dari kiri kekanan. Untuk lebih memperjelas 
masalah, segi panjang tersebut dibagi lagi menjadi beberapa lapis. 

Pada waktu t, partikel segi panjang berada dilokasi x pada permukaan. Beberapa saat 
kemudian, yaitu pada waktu t 4 At, partikel telah pindah kelokasi baru, yaitu x # V.At, 
dimana V adalah kecepatan arus. Untuk aliran potensial ini, seluruh partikel bergeser 
kekanan melalui jarak yang sama, yaitu V.At. Sekarang marilah kita periksa kasus aliran 
viskos, yang ditunjukkan pada gambar (6.8.b), (6.8.c) dan seterusnya. Lapisan paling atas, 
yang ditandai dengan arsir, dari partikel berada diluar lapisan batas atau dikawasan aliran 
potensial. Lapisan-lapisan lain dibawahnya. berada dikawasan aliran lapisan batas. Pada 
waktu t semua lapisan partikel berada dilokasi di xo dan pada waktu t # At lapisan yang 
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diarsir bergeser ke lokasi x — x9 4 Vo-At. Tetapi lapisan yang paling bawah tetap berada 
di xo karena fluida melekat pada permukaan benda. Lapisan-lapisan lainnya berada 
diantara xo dan xo # Vo-At, membentuk suatu profil yang disebut profil kecepatan. 

Seperti dapat dilihat pada gambar (6.8.b) lapisan yang di atas cenderung menyeret 
maju lapisan di bawahnya, sedangkan lapisan bawah berusaha menghambat pergerakan 
lapisan diatasnya. Sekarang marilah kita periksa keadaan dilokasi X1, disebelah hilir 
Xo,yaitu untuk xy ? xo. Partikel fluida x4 telah bergerak lebih jauh dibandingkan dengan 
partikel yang berada di xo, jadi telah mengalami hambatan viskos yang lebih besar. Lapisan 
berarsir, diluar lapisan batas, bergerak dengan kecepatan aliran potensial misalnya Vy di 
x1 dan pada waktu t # At telah bergeser ke X1 # Vg1-4t. Lapisan paling bawah tetap 
melekat di xy, sedangkan lapisan-lapisan lainnya berada dilokasi antara x1 dan xy t Vj.At. 
Karena dilokasi xy bagian partikel dalam lapisan batas telah mengalami hambatan viskos 
yang lebih besar daripada hambatan di xo, maka lebih banyak lapisan-lapisan bawah yang 
lokasinya lebih mendekati xy, yaitu profil kecepatan dibagian bawah lebih ditarik kekiri. 
Untuk partikel dilokasi xp, dimana xp » x4, lapisan-lapisan paling bawah itu lebih ditarik 
lagi kekiri, sehingga lapisan kedua dari bawah itu hampir tidak bergerak relatif terhadap 
lapisan yang melekat pada permukaan. Apabila V2? V1 ? Vp dimana V», Vj dan Vp 
adalah kecepatan di x9, x4 dan xo untuk aliran potensial diluar lapisan batas, maka tekanan 
P2 « Pj « Po karena menurut hukum Bernoulli P - 1/2pV2 - konstanta. Sebaliknya 
apabila V7 « Vj « Vg maka P7 » Pj » Po. Dibagian depan lingkaran, dari titik stagnasi 
depan sampai ketitik pundak, kecepatan aliran meningkat dengan x, yaitu V 2 Vj » Vy 
dan P2 « Pj « Po. Kondisi seperti ini disebut distribusi tekanan membantu karena selisih 
tekanan cenderung membantu mendorong lapisan partikel dalam kawasan lapisan batas 
untuk tetap bergerak maju. Sebaliknya di bagian lingkaran antara titik pundak dan stagnasi 
kecepatan fluida diperlambat kondisi tekanan yang ada adalah P» ? Pj - Po. dan disebut 
distribusi tekanan menentang. Gradien tekanan yang ada cenderung menentang atau 
mendorong lapisan partikel dalam lapisan batas kearah kebalikan arah geraknya. Untuk 
kondisi gradien tekanan membantu, partikel fluida dalam lapisan batas akan selalu tetap 
bergerak maju. Tetapi untuk kondisi gradien tekanan menentang, ada kemungkinan bahwa 
gradien tekanan yang ada itu akan mendorong partikel fluida dalam lapisan batas bergerak 
arah kebalikan dari arah semula. Semakin jauh ke hilir fluida dalam lapisan batas itu 
bergerak, semakin besarlah disipasi yang dialaminya. Hasil akhirnya adalah di lokasi x3 » 
x) beberapa lapisan ditengah itu bergerak kebalikan arah gerak fluida diluar lapisan batas 
(lihat gambar (6.6.e)).Seperti dapat dilihat lapisan kedua dari bawah bergerak kekiri, sama 
halnya dengan lapisan ketiga dan keempat. Tetapi lapisan kelima dan keenam diseret oleh 
lapisan diluar lapisan batas, dan tetap bergerak kearah kanan. Di lokasi x3 lapisan-lapisan 
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fluida yang berada dekat permukaan bergerak kekiri, sedangkan lapisan-lapisan diatasnya 
bergerak kekanan. Secara keseluruhan fluida di lapisan batas berputar membentuk suatu 
pusaran yang terpancang di kawasan sekitar x3 dan disebut gelembung pusaran (lihat 
gambar (6.8.f)). Gelembung ini jelas akan menggusur garis-garis arus secara lokal di xg. 
Garis-garis arus yang tergusur itu mungkin saja menjadi terlepas tidak mengikuti bentuk 
permukaan dan pola aliran potensial menjadi rusak sama sekali. Tetapi ada juga 
kemungkinan bahwa garis arus yang tergusur tadi kembali menempel kearah permukaan 
setelah melewati gelembung pusaran. Aliran yang menenpel kembali ini tidak begitu 
berbeda dengan aliran potensial tanpa gelembung dengan catatan bahwa lapisan batas 
disebelah hilir gelembung menjadi lebih tebal. 

Untuk kasus lain, garis arus yang tergusur menjadi terlepas sedangkan gelembung itu 
sendiri menjadi pecah. Di sebelah hilir lokasi terjadinya gelembung, suatu kawasan aliran 
alur yang acak akan terbentuk dan prediksi hasil teori aliran potensial tidak berlaku 
dikawasan alur (wake) tersebut. Dikawasan alur kecepatan aliran menjadi hampir tetap 
dengan harga yang cukup besar. Di kawasan belakang lingkaran dimana kecepatan 
potensial mengecil dan koefisien tekanan Cp seharusnya menjadi positif, ternyata Cp tetap 
berharga negatif dan besaran yang hampir tetap (lihat gambar (6.7)). Distribusi Cp untuk 
aliran yang sesungguhnya menjadi tak-simetris relatif terhadap sumbu-y. Bila distribusi Cp 
diintegralkan maka diperoleh harga koefisien gaya tahan, Cd, yang cukup besar. 

Pembahasan diatas menunjukkan bahwa gaya tahan yang dapat diamati dan diukur itu 
sesungguhnya terdiri dari 2 komponen. Komponen yang pertama adalah gaya tahan hasil 
penjumlahan tegangan geser viskos yang beraksi di seluruh permukaan, dan disebut gaya 
tahan viskos. Komponen gaya tahan yang kedua muncul karena terjadinya pelepasan aliran 
yang merusak pola aliran potensial di bagian hilir titik pelepasan. Distribusi koefisien 
tekanan dikawasan alur tersebut sangat berbeda dari distribusi hasil prediksi teori aliran 
potensial dan merusak simetri depan belakang, relatif terhadap sumbu-y, dari distribusi 
untuk aliran potensial. Oleh karena itu integral hasil penjumlahan distribusi tekanan pada 
permukaan lingkaran tidak berharga nol, dan ternyata memberikan gaya tahan yang cukup 
besar. Karena gaya tahan ini disebabkan oleh distribusi tekanan yang tak-simetris, maka 
disebut gaya tahan tekanan (pressure drag). Distribusi Cp yang tak-simetris tersebut dapat 
dilihat pada gambar (6.7). Untuk kasus aliran dengan bilangan Reynolds yang lebih besar, 
distribusi Cp ternyata lebih mendekati bentuk distribusi Cp untuk aliran potensial. Alasan 
fisis dari pengamatan ini dapat diberikan sebagai berikut. 

Untuk kasus bilangan Reynolds rendah lapisan batas dapat bersifat laminar sedangkan 
untuk bilangan Reynolds tinggi aliran lapisan batas pasti bersifat turbulen. Dalam aliran 
laminar fluida bergerak secara berlapis-lapis dan kemungkinan terjadinya pelepasan aliran 
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itu cukup besar. Dalam aliran turbulen partikel fluida bergerak acak, berputar campur- 
aduk dan mempunyai komponen kecepatan arah tegak lurus pada permukaan. Suatu 
partikel dari kawasan aliran potensial diluar lapisan batas yang mempunyai momentum 
yang besar dapat terseret masuk kedalam kawasan lapisan batas dimana partikel fluida 
mempunyai momentum yang kecil. Perpindahan momentum dari fluida Juar ke fluida 
dalam jelas membantu menambah momentum fluida dalam yang sebagai akibatnya 
terdorong bergerak lebih cepat. Dengan demikian perpindahan momentum tersebut dapat 
mengatasi hambatan tegangan geser viskos. Ini berarti bahwa proses terjadinya gelembung 
pusaran menjadi tertunda dan pelepasan aliran mungkin tidak terjadi atau bila terjadi akan 
terjadi lebih kebelakang mendekat titik stagnasi belakang. Kawasan alur yang terjadi 
menjadi semakin sempit bila titik pelepasan aliran semakin digeser kebelakang, dan 
distribusi Cp dengan demikian menjadi semakin mendekati distribusi Cp hasil prediksi 
teori aliran potensial. 

Gaya tahan tegangan geser viskos aliran turbulen memang lebih besar daripada untuk 
aliran laminar, tetapi pencegahan terjadinya pelepasan aliran jelas akan memperkecil 
komponen gaya tahan tekanan. Gaya tahan total adalah jumlah dari gaya tahan viskos dan 
gaya tahan tekanan, dan ternyata lebih besar untuk kasus aliran laminar dibandingkan 
dengan untuk aliran turbulen apabila bentuk penampang lintang silinder adalah lingkaran. 

Apabila bentuk benda itu sedemikian rupa sehingga aliran lapisan batas laminar tidak 
mudah terlepas dan menimbulkan gaya tahan tekanan yang besar, maka sebaiknya 
diusahakan supaya lapisan batas pada permukaan bersifat laminar untuk mengurangi gaya 
tahan tegangan geser viskos. Sebaiknya bila bentuk benda pasti menimbulkan pelepasan 
aliran maka sebaliknya aliran lapisan batas itu sengaja dibuat turbulen supaya pelepasan 
aliran dapat ditunda atau dicegah. Lapisan batas laminar dapat dengan sengaja dan dengan 
mudah dibuat menjadi turbulen misalnya dengan membuat permukaan menjadi sedikit lebih 
kasar atau menggunakan alat bantu pengganggu aliran seperti trip-wire dan lain 
sebagainya. 

Dalam pembahasan sebelumnya telah dinyatakan bahwa pelepasan aliran tidak akan 
terjadi dikawasan aliran dengan distribusi tekanan yang membantu. Pelapasan hanya akan 
terjadi dikawasan aliran dengan distribusi tekanan yang menentang. Pelepasan aliran juga 
akan makin lebih mudah terjadi apabila gradien tekanan yang menentang itu semakin 
curam. Sebaliknya pelepasan aliran dapat ditunda dengan membuat bentuk kurva Cp yang 
menentang itu menjadi lebih landai. 

Bentuk lingkaran memberikan distribusi Cp yang sangat membantu dari titik hidung 
ke titik pundak, tetapi sangat menentang dari titik pundak ke titik ekor. Sebaiknya bentuk 
bagian belakang benda itu dibuat lurus tidak melengkung seperti lingkaran, supaya kurva 
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distribusi Cp menentang yang terjadi juga menjadi lebih landai, tidak begitu curam. Jadi 
bentuk aerofoil yang baik adalah bentuk yang mirip lingkaran dibagian depan dan bentuk 
kurva yang lebih landai di bagian belakangnya. Bentuk aerofoil seperti itu dapat diperoleh 
dengan cara memetakan bentuk lingkaran menggunakan rumus pemetaan konformal atau 
transformasi dengan fungsi analitik variabel kompleks. Masalah pemetaan konformal ini 
akan dibahas secara lebih mendalam dalam bab berikut. 
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BAB VII 
TEORI PEMETAAN KONFORMAL 


Dalam bab sebelumnya telah dinyatakan bahwa kita ingin memetakan bentuk lingkaran 
menjadi suatu bentuk aerofoil. Pertama-tama ditentukan bahwa ada 2 bidang kompleks 
yang berbeda dan masing-masing disebut bidang-z atau bidang x-y dan bidang-C atau 
bidang € -n. Ditentukan pula bahwa pada bidang-6 terdapat lingkaran Icl 5 R yang ingin 
dipetakan menjadi suatu bentuk aerofoil di bidang-z dengan menggunakan suatu rumus 
transformasi z(C) atau € (z) yang bersifat analitik. Fungsi analitik z(C) atau € (z) itu tidak 
diketahui bentuknya dan harus dicari. 

Sebelum kita membahas masalah fungsi analitik kompleks, marilah kita tinjau ulang 
pengertian tentang fungsi analitik variabel real. Fungsi real f(x) dari variabel x adalah suatu 
hubungan matematis antara variabel x dan variabel f. Apabila rumus f(x) diberikan dan 
harga x ditentukan maka harga f dapat dihitung dengan mudah. Sebaliknya bila harga f 
diberikan maka harga x dapat dihitung. Hubungan antara harga f dan x dapat juga 
diberikan secara grafis, yaitu dengan menggunakan sistem koordinat kartesian dimana 
sumbu datarnya adalah sumbu-x, sedangkan sumbu vertikalnya adalah f. Setiap pasangan 
harga (x,f) dimana f - f(x) dapat diwakili oleh sebuah titik pada bidang kartesian tersebut. 
Apabila semua titik-titik (x,f) yang diketahui itu dihubungkan satu dengan yang lainnya, 
maka diperoleh sebuah kurva yang mewakili hubungan antara f dan x. Selanjutnya, untuk 
setiap harga x yang kita tentukan secara sembarang kita akan dapat mengetahui harga f 
yang merupakan pasangannya, dari grafik atau kurva yang diberikan. 

Sekarang marilah kita periksa apakah suatu fungsi variabel kompleks, misalnya z(6), 
dapat juga diwakili oleh suatu grafik seperti untuk fungsi real. 

Harga x dapat diwakili oleh satu garis bilangan tertentu, misalnya sumbu-x, dan 
demikian juga harga f dapat diwakili oleh suatu garis bilangan lainnya, misalnya sumbu-f, 
yang dapat ditentukan saling tegak lurus. Suatu titik kompleks G atau (E,n) tidak dapat 
diwakili oleh suatu garis bilangan tertentu, tetapi oleh bidang bilangan. Demikian juga titik 
kompleks z atau (x,y) hanya dapat diwakili oleh bidang bilangan. Pengamatan tersebut 
menunjukkan bahwa fungsi kompleks z(C) tidak dapat diwakili secara grafis oleh sebuah 
kurva seperti untuk kasus fungsi real. Tetapi ini tidak berarti bahwa fungsi kompleks z(C) 
tidak dapat diwakili secara grafis, hanya caranya saja yang berbeda. Untuk suatu harga x 
real tertentu, misalnya x - XT, suatu harga f tertentu, misalnya FT, dapat dihitung dimana 
FT - f(XT). Dalam teori fungsi variabel kompleks pernyataan bahwa untuk harga € 
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tertentu harga z — z(C) dapat dihitung itu kurang bermanfaat. Pernyataan yang lebih 
bermanfaat adalah bahwa untuk harga-harga € yang diwakili oleh titik-titik pada suatu 
kurva tertentu dibidang-C kita dapat menghitung koordinat titik-titik pada suatu kurva 
tertentu dibidang-z yang menjadi peta dari kurva dibidang-C tadi. Sebagai contoh marilah 
kita periksa fungsi kompleks berikut 


seat 11 
AD ate 7.1) 


dimana a adalah suatu konstanta real 
Untuk titik-titik yang berada di lingkaran |C| - R kita ingin mengetahui koordinat titik- 
titik pada kurva dibidang-z yang merupakan peta dari lingkaran tersebut. Apabila G-£ tin 


dan z - xtiy maka fungsi z(C) dapat diuraikan menjadi komponen real dan imajiner 
sebagai berikut 

ME mn) 5514 (2 4)) 

Ym) 11-12 417) (7.2) 
Untuk titik-titik pada lingkaran lel-R, variabel real & dan m memenuhi hubungan 
berikut 

PAPR (7.3) 
Koordinat dari titik (x,y) dibidang-z yang menjadi peta dari titik (E,n) pada lingkaran 
lel- R dapat dihitung sebagai berikut : 


aa untuk —RSEXR (7.4) 


y-n(i—(a/R)) 
Untuk kasus khusus dimana R — a, maka 


X526 

untuk —axExa (7.5) 
y-0 
Rumus diatas menunjukkan bahwa peta dari lingkaran kel 5 a adalah sepotong garis lurus, 
yaitu bagian dari sumbu-x dari x — -2a sampai ke x — 2a. Untuk kasus-kasus R - k.a, 
dimana k adalah suatu konstanta real positif, k » 1, peta dari |C| ka adalah suatu kurva 


yang koordinat titik-titiknya dapat dihitung sebagai berikut 


KK 41771 


5 | untuk —ka X€ « ka (7.6) 
y-n(k —1)/k 
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Untuk lingkaran dengan jari-jari lebih kecil daripada a, yaitu R— a/k, peta dari lcl —a/k 
adalah suatu kurva yang koordinat titik-titiknya diberikan oleh rumus berikut 


x-€(44K) 

y-n(1-K) 
Bentuk kurva yang koordinat titik-titiknya diberikan oleh rumus (7.6) dapat dipelajari 
sebagai berikut. Karena titik (£, n) berada pada lingkaran 


Pm 5 R? - (ka) (7.8) 
maka titik petanya, (x,y), berada pada kurva dengan persamaan atur berikut 


TS AAN OR NAH y 
e) “(2 p) kai 


yang dapat disederhanakan menjadi 


x 2 2 a 2 
(-) “3) 1 ON 
Titik-titik dengan koordinat yang diberikan oleh rumus (7.7) berada pada bentuk kurva 
berikut. Titik (&,n) berada pada lingkaran 
E4MP—R—(a/k) (7.10) 
jadi titik petanya, (x,y), berada pada kurva dengan persamaan atur berikut 


x “3 y 2 a 2 
(F3 Te) (2) RN 


Perhatikan bahwa rumus (7.9) itu sama dengan (7.11). Ini berarti bahwa lingkaran |C| - ka 


| untuk —(a/k)sE s (a/k) (7.7) 











dan lingkaran ll —a/k, dimana k- 1, ternyata dipetakan menjadi bentuk elips yang sama 
dibidang-z, yang persamaan aturnya diberikan oleh rumus (7.9) atau (7.11). Pembahasan 
di atas menunjukkan bahwa suatu bentuk kurva di bidang-C dapat dipetakan menjadi 
suatu bentuk kurva yang berbeda bentuk di bidang-z. Kurva di bidang-C yang bentuknya 
diketahui disebut kurva yang dipetakan, sedangkan kurva di bidang-z disebut kurva hasil 
pemetaan atau kurva peta. Rumus matematis z(C) disebut rumus pemeta atau rumus 
transformasi dari bidang-, ke bidang-z. Untuk memperoleh pengertian yang lebih 
mendalam tentang proses pemetaan, sekarang marilah kita pelajari beberapa fungsi pemeta 
sederhana. Fungsi pemeta yang paling sederhana adalah 


Z6) -6 (7.12) 
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Hasil pemetaan dengan rumus pemeta di atas adalah suatu bentuk kurva yang tepat sama 
dengan bentuk kurva yang dipetakan tanpa perubahan apapun. Karena fungsi pemeta ini 
tidak mempunyai dampak perubahan bentuk apapun maka biasanya disebut fungsi identitas 
atau fungsi tanpa dampak. Suatu fungsi pemeta sederhana yang lain adalah 

HC) -t4a (7.13) 

a—b-tic 
dimana b dan « adalah konstanta real positif. Dampak dari fungsi pemeta ini adalah 
menggeser bentuk kurva ke lokasi baru, yaitu ke arah kiri dengan jarak geser b dan ke 
bawah sebesar c, tanpa mengubah bentuknya (lihat gambar (7.1)). Fungsi pemeta (7.13) 
disebut fungsi penggeser. Fungsi pemeta sederhana lain, yang disebut fungsi pemutar, 
diberikan oleh rumus berikut 

Z6) —6.exp(-ia) (7.14) 
dimana @ adalah konstanta real positif, 

Fungsi pemutar (7.14) mempunyai dampak memutar bentuk kurva melalui sudut 
sebesar « kebalikan arah putaran jarum jam (lihat gambar (7.2)). Rumus untuk fungsi 
pengatur skala adalah 

“G)-—b.c (7.15) 
dimana b adalah konstanta real positif. Dampak dari fungsi pemeta (7.15) adalah 
mengubah skala, tanpa mengubah bentuk, kurva yang dipetakan. Apabila b » 1 maka 
fungsi (7.15) disebut fungsi pembesar skala (lihat gambar (7.3)) dan disebut fungsi 
pengecil skala apabila b « 1. 

Fungsi pemeta konstanta kompleks, yaitu 

XC) —a.t (7.16) 

a-—b-wic- Jalexp(i.argfa)) 
mempunyai dampak sebagai pengatur skala, yaitu suku modulo a atau lal, dan sebagai 
fungsi pemutar, yaitu suku argumen a atau arg(a), tetapi tidak mengubah bentuk. Fungsi 
transformasi yang diketahui itu berjumlah tak-berhingga, tetapi dampak setiap fungsi 
transformasi selalu dapat dijelaskan sebagai kombinasi dari dampak fungsi-fungsi 
transformasi sederhana yang telah dibahas sebelumnya. Fungsi transformasi yang 
diinginkan adalah fungsi yang mampu memetakan bentuk lingkaran di bidang-C menjadi 
bentuk aerofoil di bidang-z. Bentuk aerofoil adalah bentuk yang mirip lingkaran di bagian 
depan dan berbentuk kurva landai di bagian belakangnya (lihat gambar (3.1) di bab 3, 
misalnya). Rumus transformasi (7.1) telah ditunjukkan mampu memetakan bentuk 
lingkaran di bidang-4 menjadi bentuk elips, bahkan garis lurus, di bidang-z. Lingkaran 
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lcl-a dipetakan oleh rumus (7.1) menjadi garis lurus bagian dari sumbu-x dari x — - 2a 
sampai ke x — 2a. Lingkaran (c|-ka, dimana k » 1, dipetakan menjadi elips, yang 
diberikan oleh persamaan (7.9), di bidang-z. Hasil pemetaan lingkaran menjadi elips 
ditunjukkan pada gambar (7.4). Seperti dapat dilihat, untuk k — 1 atau untuk kurva a, 
lingkaran dipetakan menjadi garis. Untuk k » 1, yaitu kurva-kurva b, c dan d, lingkaran 
dipetakan menjadi elips yang menjadi semakin gemuk mendekati bentuk lingkaran untuk k 
yang makin besar. Sekarang marilah kita periksa bentuk seperti apa yang akan diperoleh 
bila rumus pemeta (7.1) diterapkan pada lingkaran yang titik pusatnya tidak berada di titik 
asal. Pertama ditentukan dulu bahwa 


8-—(k- ia 


a48 -—3(keria (7.17) 


Lingkaran yang dipetakan berjari-jari (a85) dan titik pusatnya berada di titik C, dengan 
koordinat (-8,0), atau sedikit di sebelah kiri titik asal sistem koordinat. Rumus untuk 
lingkaran ini adalah 


fe Gal (8) 4im| a45 (9.18) 
Lingkaran (7.18) bersinggungan di ujung kirinya, yaitu di titik (-(a--28),0) dengan 
lingkaran besar lel- at26. Di ujung kanan, yaitu di titik (a, 0), lingkaran (7.18) 
bersinggungan dengan lingkaran pemeta |C| — a. Ini berarti bahwa ujung kiri kurva peta 
dari lingkaran (7.18) berbentuk mirip ujung kiri bentuk elips. Di bagian kanan kurva peta 
harus menipis sehingga ujungnya mendekati bentuk garis, peta dari lingkaran kel: a. 


Seperti dapat dilihat pada gambar (7.5) hasil yang diperoleh dari pemetaan lingkaran 
(7.18) adalah sebuah aerofoil simetris. Penerapan rumus pemeta (7.2) pada titik-titik 
lingkaran memberikan koordinat titik-titik aerofoil sebagai berikut 


- (E"-5)l1 Pn AI 
“la ( kemas Eng) 


sni 5 2 - :) 
(£ -8) #(m) 
dimana €" — G—G,, yaitu 


(7.19) 
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E 533 atau in Was 
mn -n 
(FP Hn) - (a45) (7.21) 


Pusat dari lingkaran yang dipetakan tidak harus berada pada sumbu-x. Untuk kasus yang 
lebih umum, koordinat titik G, adalah (-AE,An) dan jari-jari lingkaran adalah R, 
sehingga rumus untuk lingkaran adalah 


“Ie, I-R 
ke”I-ke-6ol | (7.22) 
- —At si iAn 
Rumus di atas dapat ditulis ulang sebagai berikut 
"EA SEN 
& EtAE atau £ 8 AE 023) 
n 5-m-An atau n—-N #An 


(&-# A6) #(n- An) - R? (7.24) 
Titik (a,0) harus berada pada lingkaran, jadi jari-jari R harus memenuhi hubungan berikut 
R? — (a AE)? 4 (Any 
atau diolah lebih lanjut menjadi 


a—/R? — (An)? — AE (7.25) 


Lingkaran (7.24) bersinggungan dengan lingkaran besar di ujung kirinya, tetapi titik 
singgung berada di atas sumbu-x. Di bagian sebelah kanan lingkaran (7.24) berpotongan 
dan hampir berhimpit dengan lingkaran pemeta ll —a. Ini berarti bahwa di bagian kiri 
belahan atas bentuk aerofoil mendekati bentuk ujung kiri elips, sedangkan belahan 
bawahnya berbentuk lebih landai. Secara keseluruhan yang diperoleh adalah bentuk 
aerofoil tak-simetris (lihat gambar (7.6)). Koordinat titik-titik aerofoil diberikan oleh 
rumus berikut 








ME ym) 5 (E' sus am AH, -| 
(E” — AE) “ tAn) Fak 
Nadi Ama - a | 
yE sn (mn £ X| FE An Tar) 
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Apabila simbol G diganti dengan C” dan sebaliknya C” diganti dengan G, maka rumus- 
rumus di atas dapat diubah menjadi 








“Gm -te- ani a :) 
(£ Ke AE) ea sa An) (7.27) 
mm) - (m4 anjl 1 2 
ym) - (mt | na rar) 
atau 
e 2 Aan 
AA en (7.28) 


Go 5 —A4E tiAn 
Rumus pemeta (7.27) atau (7.28) dapat memetakan lingkaran lel- R menjadi bentuk 


aerofoil tak-simetris seperti telah dijelaskan sebelumnya. Rumus (7.28) dapat juga ditulis 
dalam bentuk berikut 


HO) XC (0) 6 Ka 


CD 561, (7.29) 
Go 5-—AE #iAn 


Rumus (7.29) dikenal sebagai transformasi Joukowski yang mampu memetakan lingkaran 
ll 5 R menjadi bentuk-bentuk aerofoil Joukowski. Bentuk aerofoil Joukowski ditentukan 
oleh 2 parameter, yaitu parameter ketebalan, (AE/R), dan parameter-parameter 


kelengkungan, (An/R). Parameter a bukan suatu parameter bebas tetapi tergantung 
pada AE dan An sebagai berikut 


a/R - /1- (An/R) — (AE/R) - (7.30) 


Suatu rumus transformasi lain yang cukup terkenal adalah rumus transformasi Karman- 
Trefftz, yaitu 
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zeka ( Uta 

z-—k.a (CG —a 

Gb, (7.31) 
Go 5—AE HiAn 


Rumus (7.31) dapat diolah menjadi bentuk z(C) atau kebaliknya, yaitu C(z), dalam 
bentuk yang lebih eksplisit sebagai berikut 


aa GA) 1G ay" 
(Gta) (CG —a)" 


z 





(7.32) 


(2-#-ka)"" 4 (z—ka)'" 
(24 ka)" — (z-— ka) 
Apabila k — 2 maka rumus (7.32) disederhanakan menjadi 





G—a (7.33) 


2 
250 4 3 
c 
jadi rumus Joukowski sesungguhnya adalah kasus khusus dari rumus Karman-Trefftz. 
Pemasukan harga k — 2 ke rumus (7.33) memberikan hasil 


€ AP ar) (9.34) 


yang merupakan kebalikan dari rumus Joukowski (7.29). Bentuk aerofoil Karman-Trefftz 
ditentukan oleh tiga parameter bentuk yaitu (AE/R), (An/R) dan k. Parameter k 
disebut parameter kelancipan ekor, karena sudut yang diapit oleh tangen kurva atas dan 
tangen kurva bawah di titik ekor, yaitu sudut t, itu ditentukan oleh harga k 

T-—(2-k).x 
Rumus Karman-Trefftz yang diterapkan untuk memetakan lingkaran kel - R, hanya dapat 


memberikan bentuk-bentuk aerofoil tertentu, bukan bentuk aerofoil yang dipilih 
sembarang. Kenyataan tersebut membatasi manfaat rumus pemeta Karman-Treffiz karena 
tidak dapat diterapkan untuk menganalisa aerofoil yang ditentukan berbentuk sembarang. 
Masalah yang dihadapi dalam analisis aerofoil berbentuk sembarang, yang ditentukan 
sebelumnya, adalah kenyataan bahwa rumus transformasinya tidak diketahui dan harus 
dicari. Salah satu cara untuk mencari rumus transformasi adalah dengan menggunakan 
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deret Laurent. Dari teori variabel kompleks diketahui bahwa di kawasan cincin (amulus) 
yang dibatasi oleh dua lingkaran konsentris di bidang-£ , yaitu di kawasan 

R atam sc Ic — Gol s Rar 
suatu fungsi analitik z(C) selalu dapat diwakili oleh deret Laurent berikut (lihat misalnya 
Wilianto (1992)) 


AG) Ba (C6) AN Ca 1G)" 35) 
dimana 
Bap aa 
“ C-6o) Ps Pp (7.36) 
Ay x9a 
NU an 


Untuk memudahkan pembahasan kita anggap bahwa titik C, berada di titik asal sistem 
koordinat, yaitu CG, - (0,0). Salah satu persyaratan yang harus dipenuhi oleh rumus 


transformasi adalah 


limit z—C (7.37) 


llbo 


Ini akan terpenuhi apabila 


B,, -0 untuk m-—2,3,... 038 
B,-1 untuk m-1 
Rumus transformasi (7.35) sekarang dapat disederhanakan menjadi 
2G) -64C, 4YC, 16" (7.39) 
m-1 
1 Ng 
Cm 5 Pa6)e” at (7.40) 
2ni 


Bila m berharga cukup besar, diperkirakan bahwa suku-suku (C,, /c") berharga 
mendekati nol dan boleh diabaikan dengan galat (error) yang kecil. Oleh karena itu untuk 
tujuan hitungan praktis, deret Laurent (7.39) dapat dipancung menjadi 


Z6) —64#C, 4G, TA oa (7.41) 
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dimana M adalah suatu bilangan bulat real positif dan bernilai cukup besar, misalnya M » 


20. Harga-harga koefisien deret Laurent, C,,, kemudian dapat dihitung dengan 


menerapkan rumus (7.40) dan menggunakan informasi tentang koordinat titik-titik pada 
aerofoil yang ingin dianalisa. 

Sekarang marilah kita anggap bahwa rumus transformasi z(C) itu sudah diketahui dan 
selanjutnya akan diterapkan untuk menganalisa masalah aliran potensial di sekeliling 
bentuk aerofoil. Fungsi potensial kompleks dari pola aliran yang ingin dianalisa harus 
memenuhi persamaan Laplace 

V?&(2)-0 (7.42) 
Dari teori mengenai fungsi analitik kompleks diketahui bahwa solusi dari persamaan 
Laplace (7.42) adalah semua fungsi analitik kompleks. Disamping itu diketahui juga bahwa 
fungsi dari fungsi analitik pasti juga bersifat analitik. Rumus transformasi C(z) bersifat 
analitik, jadi semua fungsi @(C(z)) pasti juga analitik dan merupakan solusi dari 
persamaan Laplace. Karena G(z) dapat memetakan lingkaran menjadi aerofoil atau 
sebaliknya, dan solusi untuk aliran potensial di sekeliling lingkaran adalah 


2 
00) ca si Rt (7.43) 
c 2x 
maka fungsi potensial kompleks yang dicari adalah fungsi berikut 


SN 1R 
C(z) 2x 


Apabila (&,n) adalah peta dari koordinat (x,y) yang diperoleh dari rumus transformasi 


D(2z) —6(2)- Inc(z) (7.44) 


C(z), maka harga @(z) di titik (x,y), yaitu O(x tiy), itu sama dengan harga D(C) di 
titik (E,n), yaitu O(E tim). Fungsi O(C) dan fungsi W(z) dapat diuraikan menjadi 
bentuk berikut 

DL) — DE) En) tiw(E,n) 

D(z)- O(x,y)-&(x,y) tiw(x,y) 
Karena D(z) - D(E,n) apabila (x,y) adalah pasangan titik pemetaan dari (&,n), ini 
berarti bahwa 

(5m) — &(x,y) 


Ym) —y(x,y) 
Garis arus adalah letak kedudukan titik-titik di sepanjang mana w berharga tetap. Ini 


berarti bahwa garis arus di bidang-G, misalnya w(£,n) -K, itu akan dipetakan menjadi 
garis arus w(x,y)-K di bidang-z, atau sebaliknya. Kurva isopotensial di bidang-C , 
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misalnya 6(£,n) — C, juga dipetakan menjadi kurva isopotensial 6(x,y) — C di bidang-z. 
Dalam bab sebelumnya kita telah membahas cara-cara menghitung koordinat titik-titik 
pada garis arus dan kurva isopotensial untuk aliran potensial di sekeliling lingkaran di 
bidang-C. Ini berarti bahwa garis arus dan kurva isopotensial untuk aliran di sekeliling 
bentuk aerofoil dapat juga dihitung koordinat titik-titiknya, apabila rumus pemeta z(t) 
diketahui bentuk analitiknya. Garis arus dan kurva isopotensial di bidang-C itu saling 
tegak lurus. Garis arus dan kurva isopotensial di bidang-z yang merupakan petanya juga 
saling tegak lurus. Ini berarti bahwa dua kurva di bidang-t yang saling berpotongan 
membentuk sudut @, akan dipetakan menjadi dua kurva berbentuk lain di bidang-z yang 
juga berpotongan membentuk sudut 0 yang sama. Kenyataan ini berlaku untuk kurva- 
kurva sembarang dan sudut perpotongan @ berharga sembarang, tidak terbatas pada 
pasangan kurva garis arus dan isopotensial yang saling tegak lurus saja. Bukti dari 
pernyataan tersebut tidak akan dibahas di sini, tetapi dapat dibaca di buku teks matematika 
mengenai variabel kompleks. Dapat diambil kesimpulan bahwa pemetaan dengan fungsi 
analitik kompleks itu bersifat mengawetkan sudut dan oleh karena itu disebut pemetaan 
konformal. Titik stagnasi belakang lingkaran dipetakan menjadi titik ekor aerofoil. Pada 
titik stagnasi lingkaran hanya ada satu garis singgung, sedangkan titik ekor aerofoil adalah 
titik potong dari dua kurva yang berpotongan dan mengapit sudut yang berbeda dari 180 
derajat. Ini berarti bahwa fungsi potensial kompleks, O(z), tidak bersifat mengawetkan 
sudut atau konformal pada titik ekor aerofoil, walaupun @(z) bersifat konformal untuk 
semua titik aerofoil lainnya. Kenyataan ini menunjukkan bahwa @(z) tidak bersifat 
analitik di titik ekor dan oleh karena itu turunan (d& / dz) tidak terdefinisikan di titik ekor 
aerofoil. Dari aturan turunan berrantai diketahui bahwa 

40 40 dc 40 /de Tah 

dz dc dx dt/ dt 
Dalam bab 4 telah ditunjukkan bahwa (dO/d4) adalah konjugat kompleks dari vektor 
kecepatan di bidang lingkaran atau bidang-G. Pada titik stagnasi belakang lingkaran, 
kecepatan berharga nol jadi (dO / dG) juga harus berharga nol. Supaya (d& / dz) menjadi 
tak-menentu atau tidak terdefinisikan, persamaan (7.45) menunjukkan bahwa (dz/ dt) 
harus berharga nol karena 

dD/dz-0/0 (tak-terdefinisikan) (7.46) 

Titik ekor aerofoil dimana kecepatan menjadi tak-terdefinisikan dan &(z) menjadi 

singguler, disebut sebagai titik kritis dan harus memenuhi persyaratan 
dz/ dc —0 (7.47) 
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Rumus transformasi Karman-Trefftz (7.32) dapat dicari turunannya sebagai berikut 


3 asap LEG 0 (7.48) 
a (GC ey-G-—ayT 
Koordinat titik stagnasi belakang adalah 
G5 JR? — (An) —iAn —a- AE —iAn 
dan karena €” —G #(-AE 4iAn) 
maka harga €" di titik stagnasi belakang adalah 
CG —a-tid 
Ini berarti bahwa 
C-a-0-410-0 
yang bila dimasukkan ke persamaan (7.48) memberikan hasil berikut 
dz/dt —0 
yaitu memenuhi kondisi titik kristis (7.47). Untuk menentukan kecepatan fluida di titik 
ekor aerofoil kita harus menerapkan teorema Kutta-Joukowski yang menyatakan bahwa: 








partikel fluida yang menelusuri permukaan atas aerofoil saat tiba di titik ekor harus 
mempunyai kecepatan yang sama dengan kecepatan partikel fluida yang menelusuri 
permukaan bawah aerofoil saat tiba di titik ekor. Rumus transformasi menggunakan deret 
Laurent terpancung (7.41) juga harus memenuhi kondisi titik kritis, yaitu 
LP ia, 1/3 50 (7.49) 

ac Pari 

Persyaratan (7.49) ini dapat dimanfaatkan dalam upaya menghitung harga-harga koefisien 
C,,. Dari pembahasan di atas dapat diambil kesimpulan bahwa teori pemetaan konformal 
dapat dimanfaatkan untuk mempelajari masalah aliran potensial di sekeliling bentuk 
aerofoil. Pertama-tama masalah aliran potensial di sekeliling lingkaran harus dipelajari 
terlebih dahulu. Dengan mengetahui rumus analitis dari fungsi potensial kompleks, kita 
dapat memperoleh rumus analitis untuk fungsi arus dan fungsi potensial kecepatan. 
Rumus-rumus tersebut kemudian dapat diterapkan untuk menghitung koordinat titik-titik 
pada kurva garis arus dan isopotensial. Kemudian koordinat titik-titik pada garis arus dan 
kurva isopotensial untuk aliran di sekeliling aerofoil dapat dihitung dengan menerapkan 
rumus pemeta z(C). Kecepatan aliran adalah konyugat dari turunan dari fungsi potensial 
kompleks. Apabila z(G) diketahui maka rumus untuk turunannya, yaitu dz/dt, dapat 
dicari. Distribusi kecepatan di medan aliran di sekeliling aerofoil kemudian dapat dihitung 


96 


dengan menggunakan rumus (7.45). Selanjutnya harga koefisien tekanan dapat dihitung 
dengan rumus 


C, (xy)- 1-Vx,y) (7.50) 
dimana V adalah kecepatan aliran. Sekarang ingat bahwa 


d0 
dz 


. |ao| /laz 
“ac I/ lac 


dan arah kecepatan adalah sejajar dengan tangen atau garis singgung pada permukaan 


V(xy) - (7.51) 

















aerofoil. 
Sirkulasi di sekeliling bentuk aerofoil diberikan oleh rumus berikut 





r-pv.as 

r- flag 

ii 

Tr-h a laci (7.52) 





Ini berarti bahwa sirkulasi di sekeliling lingkaran itu sama dengan sirkulasi di sekeliling 
bentuk aerofoil hasil pemetaan lingkaran tersebut dengan rumus z(C). Teori sirkulasi 
Kutta-Joukowski (lihat bab 6, persamaan (6.52)) menyatakan bahwa gaya angkat itu 
berbanding lurus dengan sirkulasi, yaitu 

L-pV.T —pV, (YV,/2R) 
Koefisien gaya angkat dirumuskan sebagai 


1 
C,- WI 20.0) 
dimana c adalah panjang tali busur. Ini berarti bahwa 
C, 5y.(4R/c) (7.53) 


dimana R adalah jari-jari lingkaran yang dipetakan, c adalah talibusur aerofoil hasil 
pemetaan dan y adalah kekuatan vorteks dalam rumus untuk @(C). 
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BAB VIII 
KARAKTERISTIK AERODINAMIK AEROFOIL 


Dalam bab sebelumnya telah ditunjukkan bahwa masalah aliran di sekeliling aerofoil 
dapat dipelajari dengan bantuan teori pemetaan konformal. Dalam bab ini pengetahuan 
tersebut akan dimanfaatkan untuk mempelajari beberapa karakteristik atau sifat 
aerodinamik bentuk aerofoil. Secara khusus kita akan mempelajari bentuk-bentuk 
aerofoil Karman-Trefftz dan sifat aerodinamiknya. Rumus pemeta Karman-Trefftz adalah 


saka Op Pele)" (8.1) 
(Ga) - (L“—a) 

(“-4-(AE —iAn) (8.2) 
dimana k, a, AE dan Am adalah konstanta real. Rumus pemetaan di atas diterapkan pada 
lingkaran 

lele --in|- R (83) 
untuk mendapatkan bentuk aerofoil Karman-Trefftz yang ingin dianalisa. Jari-jari R 
adalah fungsi dari a, AE dan An sebagai berikut 

R? — (a# A6) 4 (An)? (8.4) 
Rumus umum untuk fungsi potensial aliran angin seragam yang menyapu lingkaran 
berjari-jari R adalah 

DG) —G.e "4 Re /6 Hi(YR /21)Int (8.5) 
di mana @ adalah sudut yang diapit oleh arah angin atau arus bebas dengan arah sumbu- 
x. Turunan dari persamaan (8.5) adalah 

DG) ee Re /GHi(YR/21)/6 (8.6) 
Komponen-komponen real dan imajiner dari & adalah fungsi potensial kecepatan, &, dan 
fungsi arus, W. Rumus untuk & dan w dalam koordinat kartesian adalah 








sentana gtaj Ba) on 
VE (neon Aaj Ru 69 
E4tn 21 


Rumus rumus di atas untuk koordinat polar adalah 


&(r,0) — ( 4# 3 Jeese @«)-— In 
r 21 





0 (8.9) 


98 





vso-|r 8 Jimce Oa (8.10) 
r 2x 


Turunan fungsi & adalah konjugat kompleks dari vektor kecepatan, jadi rumus untuk 
komponen kecepatan dapat diperoleh dari persamaan (8.6). Dalam koordinat kartesian, 
komponen-komponen kecepatan adalah 











R' (8-1) R'En Rn 
En) Il 2 ina 8.11 
2 | Ta Na aa Ne Ta agar 14 
R1). R'En & £ 
sm) 5I1- 2 —— 8.12 
TAN | (Snp PN amp Em ai 
Dalam koordinat polar, komponen-komponen kecepatan diberikan sebagai berikut 
R' 
V, (1,0) | 1- — Icos(0-—“) (8.13) 
r 
2 
vatn0) (1s Janto-a -e (8.14) 
Pada permukaan lingkaran, yaitu untuk r — R 
V, (R,0) -—0 (8.15) 
Vo (R,0) - -2sin(O —a)- 1 (8.16) 


Seperti dapat dilihat, kecepatan partikel fluida itu arahnya sejajar dengan tangen atau 
garis singgung pada permukaan lingkaran. Ini sesuai dengan kondisi batas yang menuntut 
bahwa fluida tidak boleh menembus keluar masuk permukaan silinder bundar. 

Letak titik stagnasi, dimana kecepatan berharga nol, dapat diperoleh dari 
persamaan (8.16) sebagai berikut 

sin(0,, —«) ——y/4x 


Bila dirumuskan bahwa 


B-sin“ (1/41) (8.17) 
maka koordinat titik stagnasi depan adalah (R,0,,, ) dimana 

Oa 5atT-4P (8.18) 
sedangkan koordinat titik stagnasi belakang adalah (R,0,,,) dimana 

Os —a—P (8.19) 


Rumus untuk garis arus diberikan oleh rumus (8.8) atau rumus (8.10), yang dapat diolah 
lebih lanjut sebagai berikut. Untuk fungsi arus bernilai K (konstanta), rumus (8.10) dapat 
ditulis ulang menjadi 
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sin(@—a) — (K/R-(y/217)Inr)/(T/R-R/r) (8.20) 
dimana &, K, R dan y diketahui nilai-nilainya. Apabila harga diberikan, dengan catatan 
bahwa r» R, maka harga @ dapat dihitung dan koordinat titik (r,@) yang berada pada 
garis arus W — K diperoleh 

Pola aliran disekeliling lingkaran untuk kasus « - O ditunjukkan pada gambar (8.1) 
dan untuk « #0 diberikan pada gambar (8.2). Lingkaran |(c|-R dapat dipetakan 


menjadi bentuk aerofoil oleh rumus pemetaan Karman-Trefftz (8.1) seperti telah 
dijelaskan dibab sebelumnya. Perlu diingat bahwa titik ekor aerofoil adalah peta dari titik 
stagnasi belakang lingkaran. Selanjutnya dapat diperhatikan juga bahwa rumus pemeta 
Karman-Trefftz menuntut agar koordinat titik stagnasi belakang lingkaran harus 
diberikan dalam sistem koordinat Kartesian oleh rumus berikut 


(65m). — (at AE,—An) (8.21) 
atau dalam sistem koordinat polar oleh rumus 
(R,0,) - (R,-sin” (An/R)) (8.22) 


Perbandingan antara rumus (8.22) dan (8.19) menunjukkan bahwa 
sinf(a — B)—- —An/R 
Bila dirumuskan bahwa 


sin2 — An/R (8.23) 
maka rumus sebelumnya dapat ditulis menjadi 

«-B41—0 (8.24) 
Rumus (8.24) dapat ditulis ulang sebagai 

P-atA (8.25) 


Dari bab 6 diketahui bahwa gaya angkat pada silinder bundar diberikan oleh rumus 
berikut 


C, -2y (8.26) 
Selanjutnya dari rumus (8.17) diketahui bahwa 

y —4nxsinp (8.27) 
Gabungan rumus (8.25), (8.26) dan (8.27) memberikan hasil akhir sebagai berikut 

C, —- 8nsinf(a#A) (8.28) 


Dalam bab 7 telah ditunjukkan bahwa (An/R) adalah parameter kelengkungan aerofoil. 
Apabila An — 0 maka lingkaran dipetakan menjadi aerofoil simetris dan rumus (8.28) 
disederhanakan menjadi 

C, — 8xsina (8.29) 
Pengamatan diatas menunjukkan bahwa gaya angkat untuk aerofoil simetris itu 
berbanding lurus dengan sina. Untuk @ bernilai kecil, dalam satuan radian, diketahui 
bahwa 

sina sa (8.30) 


jadi gaya angkat kira-kira berbanding lurus dengan @ untuk @ bernilai cukup kecil. 
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Pembahasan bab 7 juga menunjukkan bahwa koefisien gaya angkat untuk aerofoil hasil 
pemetaan lingkaran |C| - R diberikan oleh rumus berikut 

C, —y(4R /e) (8.31) 
dimana c adalah talibusur aerofoil. 
Ini berarti bahwa untuk aerofoil simetris 


C, -16x(R /c)sina (8.32) 
dan khususnya bila « bernilai cukup kecil maka 
C, s167(R/c)a (8.33) 


Untuk Am - 0 atau untuk aerofoil simetris, rumus pemeta Karman-Trefftz (8.1) dapat 
disederhanakan menjadi bentuk berikut 


eka GA ta)" 4 (C— AE —a)" 
(G- AE ta)" — (C- AE —a)" 


atau karena R-a-AE, maka 








ka GER- 246)" 4 (C-R)" 
(GR - 246)" — (G-R)" 
Titik ekor aerofoil adalah peta dari titik stagnasi belakang lingkaran, yaitu titik (R,0), jadi 
koordinat titik ekor adalah 


ai 2AE,0)" 4 (R-R,0)" 
(2R — 2AE,0)" — (R- R,0)" 
yaitu titik dengan koordinat (ka,0). 
Titik hidung aerofoil adalah peta dari titik lingkaran (-R,0), jadi koordinat hidung 
aerofoil adalah 








3 ka 246,0)" 4 (-2R,0)" 
(-24E, 0)" — (-2R, 0)" 
atau karena —1 — exp(—it) maka 


z 


ka (zat) exp(—ikt) t lari" exp(-ikt) 
ec 
Izat|" exp(-ikx)— (2R|" exp(—ikn) 





(ME /R)" —1” 
Panjang talibusur adalah jarak dari titik hidung ke titik ekor, jadi panjang talibusur c 
adalah 


a-| aa BR, #1 ) 
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(ME/R) 41 


c— ka 
(ME/R)F —1 





yaitu 
c- 2ka /JL-— (AE /R)" 
c-2k(R- AE)/(I- (AE/R)") 
Hasil diatas dapat ditulis ulang sebagai 


1-— (AE/R)" 


R/c- (8.34) 
2k(1- AE/R) 
Gabungan rumus (8.32) dan (8.34) memberikan hasil berikut 
e k 
Cc, 2 SELAMEARY na: (8.35) 


k 1-AE/R 

Apabila k — 2, yaitu untuk aerofoil Joukowski simetris, rumus (8.35) dapat 
disederhanakan menjadi 

C, 5-4n(14#- AE /R)sina (8.36) 
Seperti dapat dilihat harga koefisien gaya angkat ternyata meningkat dengan harga 
AE/R, yaitu dengan parameter ketebalan. Karena tebal aerofoil meningkat secara 
monoton dengan harga (A£ /R), maka dapat diambil kesimpulan bahwa koefisien gaya 
angkat meningkat dengan tebal aerofoil. 
Untuk A£ — 0, rumus (8.35) dapat disederhanakan menjadi 


C,- “Asing (8.37) 


Ini berarti bahwa koefisien gaya angkat mengecil dengan semakin membesarnya harga k. 
Karena sudut kelancipan ekor aerofoil mengecil dengan semakin membesarnya harga k, 
maka koefisien gaya angkat akan membesar dengan semakin membesarnya sudut 
kelancipan ekor. Walaupun analisis diatas dilakukan untuk aerofoil Karman-Treffiz, 
tetapi pada dasarnya berlaku juga secara umum. 

Kesimpulan umum yang dapat diambil adalah bahwa koefisien gaya angkat menjadi 
semakin besar untuk aerofoil yang lebih tebal dan sudut kelancipan ekornya lebih besar. 
Sekarang perhatikan bahwa rumus (8.28) untuk sudut serang berharga nol, yaitu « — 0, 
dapat disederhanakan menjadi 


C, — 8xsin/ 
yang dapat digabung dengan rumus (8.23) menjadi 
C, -87(An/R) (8.38) 


Ini berarti bahwa koefisien gaya angkat berharga lebih besar untuk harga parameter 
kelengkungan, An/R, yang membesar. Karena kelengkungan aerofoil meningkat 
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dengan meningkatnya harga An/R, maka koefisien gaya angkat aerofoil pada sudut 
serang tertentu akan meningkat untuk aerofoil dengan kelengkungan yang lebih besar 
Untuk mempelajari pola aliran di sekelililng suatu bentuk aerofoil Karman-Trefftz, 
kita harus melakukan langkah-langkah berikut. Pertama-tama kita harus menentukan 
harga-harga parameter bentuk, yaitu parameter ketebalan, AE/R, parameter 
kelengkungan, An/R, dan parameter kelancipan ekor, k. Harga-harga tersebut 
menentukan bentuk aerofoil yang ingin diselidiki. Selanjutnya kita tentukan harga sudut 
serang &, yaitu sudut antara arah angin atau arus bebas dengan arah talibusur aerofoil. 
Harga sirkulasi, y, kemudian dihitung dengan rumus (8.17) yang dapat ditulis ulang 


menjadi 

y 54mxsinp (8.39) 
dimana harga B dihitung dengan rumus (8.25), yaitu 

P-atA 


dan A dihitung dengan rumus (8.23), yaitu 
Asin" (An/R) 

Pola aliran di sekeliling lingkaran, dengan harga « yang diberikan dan harga y yang 
dihitung dengan rumus (8.39), kemudian dihitung dengan menggunakan rumus (8.8) atau 
(8.10). Koordinat titik-titik stagnasi depan dan belakang lingkaran kemudian dihitung 
menggunakan rumus (8.17), (8.18), dan (8.19). Koordinat titik-titik lingkaran kemudian 
dipetakan menjadi koordinat titik-titik permukaan aerofoil. Titik stagnasi depan lingkaran 
dipetakan menjadi titik stagnasi depan aerofoil, sedangkan titik stagnasi belakang 
lingkaran dipetakan menjadi titik ekor aerofoil. Titik-titik di sepanjang garis arus pada 
pola aliran di sekeliling lingkaran kemudian dipetakan menjadi titik-titik di sepanjang 
garis arus pada pola aliran di sekeliling aerofoil. Kecepatan aliran di titik (x,y) pada 
permukaan aerofoil yang merupakan peta dari titik (&,n) pada lingkaran dapat dihitung 
sebagai berikut. Pertama-tama kecepatan di titik lingkaran (£,n) dihitung menggunakan 


rumus berikut 


V(£,n) —- V(R,0) - 2sin(0 —a)-#-y/21x (8.40) 
O-—tan" (n/E) 
Kecepatan di titik aerofoil (x,y) kemudian dihitung dengan rumus berikut 
VGsy) — V(£,m) /|dz/ dc| (8.41) 


x Xx. k-1 
IE oyay AE HAN Era 
da, NC “ta)" —(C“-—a)I 
C#- (E-—AE)ti(ntAn) 
Setelah harga V(x,y) diperoleh, kemudian harga koefisien tekanan di titik (x, y), yaitu 
C, (x,y), dihitung sebagai berikut 
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C, (x,y) 51-V (x,y) (8.42) 
Pola-pola aliran di sekeliling bentuk aerofoil lengkung untuk kasus sudut serang « — 0 
ditunjukkan pada gambar (8.3), sedangkan untuk « # 0 ditunjukkan pada gambar (8.4). 
Dapat diamati bahwa talibusur aerofoil hasil pemetaan itu selalu sejajar dengan sumbu-x, 
sedangkan arah arus bebas membentuk sudut d« dengan arah sumbu-x. Profil distribusi 
koefisien tekanan, C,, sebagai fungsi ordinat aerofoil, x, untuk kasus sudut serang d« — 0 
ditunjukkan pada gambar (8.5) dan untuk « #0 ditunjukkan pada gambar (8.6). 
Perhatikan bahwa sumbu vertikal untuk C, itu telah diputar balik dari biasanya, yaitu C, 
bernilai negatif untuk bagian di atas sumbu-x dan positif di bawah sumbu-x. Cara 
penyajian data C, seperti itu didasarkan pada kenyataan bahwa tekanan negatif 
cenderung menyedot benda ke arah kawasan dimana tekanan bernilai negatif. Distribusi 
C, dapat digunakan untuk menghitung koefisien gaya angkat dan arah gaya angkat 
positif adalah ke atas, tegak lurus pada arah arus bebas. Cara penyajian data C, seperti 
dijelaskan di atas memberikan kesan visual mengenai arah gaya angkat. Disamping itu 
harga koefisien gaya angkat dapat ditaksir dengan cepat secara visual, yaitu sebagai luas 
kawasan yang dikurung oleh sirkit atau kurva tertutup profil distribusi C,: 

Perhatikan juga bahwa data koordinat titik-titik pada aerofoil itu telah 
dinormalisasikan, yaitu dibagi oleh panjang talibusur aerofoil. Dengan cara ini semua 
aerofoil selalu disajikan sebagai mempunyai talibusur sebesar satu satuan panjang. 
Selanjutnya letak aerofoil selalu digeser, diatur sedemikian rupa sehingga titik ekor 
aerofoil mempunyai koordinat (1,0). Pengaruh ketebalan aerofoil pada bentuk profil C, 
dipelajari dengan memeriksa aerofoil Karman-Trefftz dengan harga-harga An,a dan k 
tertentu dan memvariasikan harga AE. Hasil yang diperoleh ditunjukkan pada gambar 
(8.7) untuk kasus aerofoil simetris, yaitu untuk An — 0 dan sudut serang &« — 0. Untuk 


kasus tersebut gaya angkat bernilai nol dan ini tercermin pada kenyataan bahwa profil C, 
untuk belahan atas aerofoil itu tepat sama dan berhimpit dengan profil C,, untuk belahan 


bawah aerofoil. 
Untuk harga A& yang membesar, atau untuk aerofoil yang semakin tebal, harga 


mutlak C, maksimal juga semakin besar. Letak titik C, maksimal semakin mundur ke 
belakang, ke arah ekor dan profil C, pada bagian belakang C, maksimal menjadi 
semakin curam. Karena distribusi tekanan di bagian belakang adalah distribusi tekanan 
yang menentang, maka profil C, yang makin curam itu kurang baik dipandang dari segi 
kemungkinan terjadinya pelepasan aliran sebagai akibat adanya lapisan batas dimana 


viskositas berdampak besar. Pengaruh kelengkungan aerofoil pada profil C, dapat dilihat 
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pada gambar (8.8). Harga-harga parameter pemetaan Karman-Trefftz, yaitu AE, a dan k, 
ditentukan dan @ dipilih bernilai nol sedangkan harga An divariasikan. Karena An # 0 


maka aerofoil hasil pemetaan berbentuk tak-simetris dan profil C, untuk belahan atas 


aerofoil berbeda dari profil C, untuk belahan bawahnya. Untuk kasus ini gaya angkat 
tidak berharga nol dan besarannya diberikan secara visual oleh luas kawasan yang 
dikurung oleh sirkit atau kurva tertutup profil C,, untuk belahan atas yang menyambung 


dengan profil untuk belahan bawah. Seperti dapat dilihat, luas tersebut semakin besar 
untuk harga An yang semakin besar, yaitu untuk bentuk aerofoil yang semakin lengkung. 


Harga C, mutlak maksimal untuk belahan atas meningkat dengan semakin membesarnya 
harga Am. Letak titik C, maksimal mundur ke belakang dengan meningkatnya harga 
Am, tetapi pada dasarnya profil C, belahan atas itu tidak terlalu curam atau cukup 


landai. Sebaliknya profil C, untuk belahan bawah menjadi semakin curam dengan 
meningkatnya harga An, sehingga kemungkinan terjadinya pelepasan aliran di belahan 
bawah itu cukup besar untuk An yang besar. 

Parameter bentuk terakhir yang perlu dipelajari adalah parameter kelancipan ekor, k. 
Profil distribusi C,, untuk beberapa aerofoil dengan harga AE, Am dan a tertentu, dan k 
divariasikan ditunjukkan pada gambar (8.9). Dapat diamati bahwa untuk variasi harga k 
yang relatif besar, profil C, yang diperoleh ternyata hanya sedikit berbeda. Dengan 
semakin mengecilnya nilai k, atau sudut kelancipan ekor semakin membesar, harga 


mutlak C, cenderung meningkat baik untuk belahan bawah maupun belahan atas 
aerofoil. Letak titik harga mutlak C, maksimal cenderung mundur ke belakang. Di 
bagian belakang aerofoil, profil C, cenderung menjadi semakin curam. Luas kawasan 


yang dikurung oleh sirkit profil C, cenderung bertambah luas, yaitu gaya angkat 
meningkat dengan semakin besarnya sudut kelancipan ekor. Pengamatan ini sesuai 
dengan prediksi harga C, dengan menggunakan rumus (8.35) atau (8.37). Dari 


pengamatan-pengamatan tersebut di atas kita dapat mengambil kesimpulan sebagai 
berikut. Untuk sudut serang « —0 bentuk aerofoil yang paling baik, dipandang dari segi 
gaya angkat yang maksimal dan menghindari terjadinya pelepasan aliran, adalah aerofoil 
dengan ketebalan sedang, cukup lengkung untuk belahan atasnya tetapi landai pada 
belahan bawahnya dan sudut kelancipan ekor yang sedang atau kecil. Walaupun 
kesimpulan di atas diambil berdasarkan pengamatan pada aerofoil Karman-Trefftz, tetapi 
garis besarnya berlaku secara umum. 


Sekarang marilah kita pelajari pengaruh sudut serang « pada profil distribusi C, 


untuk sebuah aerofoil Karman-Trefftz tak-simetris tertentu. Bentuk profil C, tersebut 
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dapat dilihat pada gambar (8.10). Pada belahan bawah aerofoil, peningkatan sudut serang 
cenderung mengubah bagian depan bentuk profil C, yaitu menariknya ke bawah, 
sedangkan bagian belakangnya cenderung terangkat ke atas walaupun untuk bagian 


belakang ini perubahan bentuk profil C, sangat kecil. Perubahan yang jauh lebih besar 
terjadi pada belahan atas aerofoil. Harga C, mutlak maksimal cenderung meningkat 


pesat untuk nilai « yang cukup besar dan letak titik C, maksimal tersebut maju ke depan 
untuk « yang meningkat. Ini berarti bahwa peningkatan harga « diikuti oleh 
meningkatnya kecuraman profil C,, yang tentu saja diiringi meningkatnya kemungkinan 
terjadinya pelepasan aliran. Lebih parah lagi letak dari titik pelepasan aliran ini bergerak 
maju ke depan untuk & yang membesar, jadi hampir seluruh permukaan atas aerofoil 
akan berada di kawasan alur yang tidak diinginkan. Memang gaya angkat aerofoil 
cenderung meningkat dengan meningkatnya harga @ , jadi gaya angkat yang besar akan 
diperoleh untuk « yang besar. Tetapi untuk « yang cukup besar, pelepasan aliran pasti 
terjadi di permukaan atas dan gaya angkat menurun secara mendadak. Fenomena ini 
disebut stall atau aliran mogok dimana gaya angkat menurun sedangkan gaya tahan 
meningkat secara mendadak. Masalah ini tidak dapat diramalkan dengan menggunakan 
teori aliran potensial, tetapi harus dipelajari dengan menggunakan model aliran viskos 
yang terlepas dari permukaan aerofoil dan membentuk aliran pusaran atau alur. Dari teori 
lapisan batas viskos, diketahui bahwa lapisan batas dimana viskositas tidak boleh 
diabaikan itu dalam kenyataannya sangat tipis, kira-kira kurang daripada seperseribu 
panjang talibusur aerofoil. Di luar lapisan batas, pola aliran dapat diramalkan dengan baik 
oleh hasil-hasil perhitungan teori aliran potensial. Oleh karena itu untuk aliran yang tidak 
terlepas, teori aliran potensial yang diberi koreksi atau perbaikan untuk dampak lapisan 
batas, dapat diterapkan untuk menghitung gaya angkat dan gaya tahan aerofoil dengan 
hasil yang cukup memuaskan. Untuk dapat menerapkan koreksi lapisan batas, kita 
membutuhkan informasi mengenai profil disribusi kecepatan aliran di sepanjang 
permukaan aerofoil dari titik stagnasi depan sampai ke titik ekor. Dsitribusi kecepatan 
aliran di sepanjang permukaan silinder bundar (lingkaran) diberikan oleh rumus (8.16) 
dan ditayangkan pada gambar (8.11) untuk y — O dan gambar (8.12) untuk y #0 . 

Distribusi kecepatan aliran di sepanjang permukaan aerofoil diberikan oleh rumus 
(8.41) dan ditayangkan pada gambar (8.13) dan (8.14) untuk aerofoil simetris dan pada 
gambar (8.15) dan (8.16) untuk aerofoil tak-simetris. 
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BAB IX 


PENUTUP 


Dalam buku ini kita telah membahas konsep-konsep dasar dinamika fluida. Walaupun 
sesungguhnya fluida itu terdiri dari partikel-partikel kecil seperti atom dan molekul, tetapi 
pembahasan dinamika fluida hanya dapat dilakukan bila kita menganggap fluida sebagai 
materi kontinyu. Konsep tentang anggapan materi kontinyu dan partikel fluida telah 
dijelaskan secara rinci. Pergerakan fluida materi kontinyu dapat dipelajari dengan 
menerapkan beberapa teori dasar Fisika, yaitu hukum-hukum kelanggengan massa, 
momentum dan energi. Pergerakan fluida dapat dipelajari dengan menggunakan 2 
pendekatan, yaitu pendekatan Euler dan pendekatan Lagrange. Dalam buku ini hanya 
pendekatan Euler yang dibahas. 

Persamaan matematis yang mengatur pergerakan fluida dalam bentuk paling lengkap 
yang dibahas dalam buku ini adalah persamaan Navier-Stokes. Apabila viskositas fluida 
dianggap bernilai nol maka persamaan Navier-Stokes dapat disederhanakan menjadi 
persamaan Euler. Dengan menganggap bahwa densitas fluida bernilai tetap dan aliran 
fluida bersifat tak-rotasional (tanpa pusaran) maka persamaan Euler selanjutnya dapat 
disederhanakan menjadi persamaan Laplace. Dalam kenyataan viskositas fluida bernilai 
kecil tetapi bukan nol dan untuk banyak kasus dampak viskositas tidak boleh diabaikan. 
Oleh karena itu pembahasan singkat mengenai pengaruh viskositas telah diberikan dalam 
buku ini. Namun demikian karena masalah aliran fluida viskos itu sangat rumit maka dalam 
buku ini pembahasan lebih banyak dipusatkan pada aliran fluida potensial yang tak-viskos. 
Persamaan atur untuk aliran potensial adalah persamaan Laplace, jadi pokok bahasan buku 
ini adalah teori dasar dan metoda penyelesaian persamaan Laplace. Persamaan ini dapat 
mudah dimengerti dengan menggunakan pendekatan metoda variabel kompleks. Oleh 
karena itu dalam buku ini teori variabel kompleks dibahas secara cukup mendalam untuk 
memperoleh gambaran yang jelas mengenai solusi matematis dari persamaan Laplace dan 
sifat-sifat dari solusi-solusi tersebut. Pendekatan matematis memang diperlukan untuk 
pengolahan masalah lebih lanjut, tetapi dinamika fluida tidak sama dengan matematika. 
Dalam mempelajari dinamika fluida, yang diinginkan adalah pengertian fisis mengenai 
pergerakan partikel fluida dan interaksinya dengan permukaan benda padat yang 
membatasi ruang gerak fluida tersebut. Untuk memperoleh pengertian fisis tersebut kita 
telah membahas secara cukup mendalam mengenai beberapa konsep dasar dalam 
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pergerakan fluida, seperti garis arus, fungsi arus, kurva isopotensial, fungsi potensial, 
vortisitas, sirkulasi, koefisien tekanan, koefisien gaya angkat dan lain sebagainya. Solusi- 
solusi elementer persamaan Laplace seperti sumber/serap, vorteks dan doblet yang 
diketahui bentuk persamaan matematiknya telah dicoba diberi penjelasan fisis yang 
diharapkan dapat memberikan pengertian yang lebih mendalam. Sifat linearitas persamaan 
Laplace dimanfaatkan untuk memperoleh solusi-solusi lain yang lebih rumit sebagai hasil 
kombinasi linear solusi-solusi elementer yang sederhana tadi. Solusi yang lebih rumit ini 
dipaksa memenuhi syarat batas yaitu mewakili masalah aliran yang lebih nyata dan 
bermakna, seperti aliran di sekeliling silinder bundar. Aliran ini dapat diperoleh dengan 
mudah sebagai kombinasi linear aliran angin seragam, doblet dan vorteks. Walaupun aliran 
di sekeliling lingkaran (penampang silinder bundar) itu cukup realistis dipandang dari segi 
fisis, tetapi tidak begitu bermanfaat dari segi rekayasa. Masalah yang sesungguhnya ingin 
dipelajari adalah pola aliran di sekeliling bentuk aerofoil, yang kemudian diharapkan dapat 
memberikan pengertian fisis mengenai karakteristik aerodinamik bentuk-bentuk aerofoil. 
Untuk mencapai tujuan tersebut kita telah membahas masalah teori pemetaan 
konformal dengan variabel kompleks, yang memberikan kemampuan untuk memetakan 
pola aliran di sekeliling lingkaran menjadi pola aliran di sekeliling aerofoil. Dalam buku ini 
rumus pemeta yang dibahas secara rinci adalah rumus yang sederhana yaitu transformasi 
Joukowski dan transformasi Karman-Trefftz. Manfaat kedua rumus pemeta tersebut 
memang agak terbatas karena tidak dapat diterapkan pada bentuk aerofoil sembarang. 
Namun demikian rumus-rumus tersebut dapat digunakan sebagai wahana untuk 
menjelaskan bagaimana informasi tentang pola aliran, distribusi kecepatan dan koefisien 
tekanan untuk aliran di sekeliling lingkaran dapat diolah lebih lanjut menjadi informasi 
tentang aliran di sekeliling bentuk aerofoil. Kedua rumus transformasi tersebut juga 
bermanfaat karena mampu mengidentifikasikan 3 parameter bentuk aerofoil yaitu 
parameter ketebalan, kelengkungan dan kelancipan ekor aerofoil. Analisis tentang 
karakteristik aerodinamik bentuk aerofoil Joukowski dan Karman-Treffiz ternyata juga 
cukup bermanfaat untuk mendapatkan pengertian mengenai bagaimana distribusi koefisien 
tekanan dan gaya angkat sebuah aerofoil di pengaruhi oleh ketebalan, kelengkungan dan 
kelancipan ekornya. Analisis tersebut juga mampu memberikan informasi tentang 
bagaimana koefisien gaya angkat sebuah bentuk aerofoil bervariasi dengan sudut 
serangnya. Pada dasarnya bahan-bahan yang dibahas secara rinci dalam buku ini cukup 
untuk memberikan pengertian yang mendalam tentang pergerakan fluida tak-mampat, tak- 
viskos, tak-rotasional, tunak dan bersifat 2 dimensional. Analisis tentang karakteristik 
aerodinamik aerofoil dalam aliran potensial tersebut juga dapat dilakukan secara tuntas 
walaupun terbatas pada bentuk-bentuk aerofoil Joukowski dan Karman-Trefftz. Untuk 
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menganalisa bentuk aerofoil sembarang, kita membutuhkan suatu rumus transformasi yang 
bersifat lebih umum seperti transformasi Deret Laurent. Konsep dasar transformasi dengan 
deret Laurent juga telah diberikan, walaupun tidak secara terperinci. Ruang lingkup buku 
ini memang dibatasi hanya sebagai pengantar atau pengenalan awal. Masalah analisis dan 
perancangan aerofoil dengan banyak komponen berbentuk sembarang memang belum 
dibahas secara rinci di buku ini, tetapi menjadi pokok bahasan dalam modul kedua bahan 
kursus ini. Dalam modul kedua tersebut masalah transformasi dengan deret Laurent akan 
dibahas secara tuntas. Disamping itu konsep singgularitas kontinyu yang mendasari 
metoda panel dan teori tentang metoda pancl untuk menganalisa dan merancang bentuk- 
bentuk acrofoil juga akan dibahas secara rinci dalam modul kedua tersebut, 








Gbr. 6.3 Distribusi Cp sepanjang sumbu-x 











Gbr. 6.5 Aliran disekeliling lingkaran dengan 


sirkulasi y - 7/2 











Gbr. 6.6 Distribusi kecepatan dan Cp disekeliling 
lingkaran dengan sirkulasi y - 1/2 





Gbr. 6.7 Perbandingan distribusi Cp hasil teori dan eksperimen 
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Gbr. 6.8 Dampak viskositas dan gelembung resirkulasi 
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Gbr. 7.1 Dampak fungsi penggeser 
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Gbr. 7.2 Dampak fungsi pemutar 
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Gbr. 7.3 Dampak fungsi pengatur skala 
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Gbr. 7.4 Transformasi Joukowski 











Gbr. 7.5 Aerotoil simetris Joukowski 











Gbr. 7.6 Aerofoil lengkung Joukowski 





Gbr. 8.2 Pola aliran disekeliling lingkaran. « # 0 








Gbr. 8.3 Pola aliran disekeliling aerofoil lengkung 
hasil pemetaan pola aliran pada gambar 8.1 
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Gbr. 8.4 Pola aliran disekeliling aerofoil simetris, sudut 
serang &, hasil pemetaan pola aliran pada gambar 8.2 
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Gbr. 8.5 Distribusi Cp (sudut serang « — 0) 
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Gbr. 8.6 Distribusi Cp (sudut serang c0) 
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Gbr. 8.7 Distribusi Cp, aerofoil simetris, « - 0 
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Gbr. 8.8 Distribusi Cp, aerofoil lengkung. « — 0 
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Gbr. 8.9 Distribusi Cp, aerofoil lengkung, « — 0 
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Gbr. 8.10 Distribusi Cp, aerofoil simetris 
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Gbr. 8.11 Distribusi kecepatan (y — 0) 





Gbr. 8.12 Distribusi kecepatan (y #0) 
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Gbr. 8.13 Distribusi kecepatan sepanjang 
permukaan aerofoil simetris, « — 0 
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Gbr. 8.15 Aerofoil lengkung, « -0 
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Gbr. 8.16 Aerofoil lengkung, « —5 
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Gbr. 2.1 Pergerakan dan kesetimbangan gaya 
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Gbr. 3.1 Kesetimbangan gaya pada pesawat dan 
pola aliran disekeliling aerofoil atau sayap 








Gbr. 3.2 Aliran dalam pipa planar, 2-dimensional 








Gbr. 3.3 Aliran dalam pipa 1-dimensional semu 
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Gbr. 3.4 Fluks (laju aliran) massa antara 2 titik 
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Gbr. 3.5 Fungsi arus, garis arus dan kecepatan aliran 
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Gbr. 3.6 Kinematika pergerakan partikel fluida 





Gbr. 3.7 Gerakan berputar (rotasi) benda kokoh 
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Gbr. 3.8 Tornado dan distribusi kecepatan untuk aliran vorteks 





Gbr. 3.9 Sirkulasi disekeliling sirkit C 
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Gbr. 3.10 Sirkit C bentuk sembarang dibagi menjadi 
sejumlah besar elemen segi panjang 





Gbr. 3.11 Sirkulasi pada sirkit ABCTA 
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Gbr. 4.1 Bidang kompleks dalam sistem koordinat kartesian 
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Gbr. 4.2 Diagram argand 
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Gbr. 5.1 Aliran angin seragam 
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Gbr. 5.2 Aliran sumber/serap 
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Gbr. 5.3 Penjelasan fisis aliran sumber/serap 
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Gbr. 5.4 Garis sumber 2-D dan titik-sumber 3-D 
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Gbr. 5.5 Aliran vorteks 


lubang pada 
dasar bejana 


Gbr. 5.6 Vorteks dalam bejana berlubang didasarnya dan berisi air 
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Gbr. 5.7 Sirkit-sirkit tertutup yang tidak mengurung kawasan inti vorteks 





Gbr. 5.8 Aliran doblet 







































































Gbr. 6.1 Aliran disekeliling bentuk lingkaran 


